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In these lecture notes we review some recent attempts at searching for non-Fermi liquids and 
novel quantum phase transitions in holographic systems using gauge/gravity duality. We do this 
by studying the simplest finite density system arising from the duality, obtained by turning on a 
nonzero chemical potential for a U(l) global symmetry of a CFT, and described on the gravity side 
by a charged black hole. We address the following questions of such a finite density system: 

1. Does the system have a Fermi surface? 

What are the properties of low energy excitations near the Fermi surface? 

2. Does the system have an instability to condensation of scalar operators? 
What is the critical behavior near the corresponding quantum critical point? 

We find interesting parallels with those of high T c cuprates and heavy electron systems. Playing a 
crucial role in our discussion is a universal intermediate-energy phase, called a "semi-local quantum 
liquid," which underlies the non-Fermi liquid and novel quantum critical behavior of a system. It 
also provides a novel mechanism for the emergence of lower energy states such as a Fermi liquid or 
a superconductor. 
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I. 



INTRODUCTION AND MOTIVATION 



Understanding phases of matter for which there is no quasiparticle description presents some of the most challenging 
problems in physics. Examples include the quark-gluon plasma (QGP) created in heavy ion collisions at RHIC during 
the past decade and now at LHC and ultracold atomic systems at unitarity. In condensed matter physics prominent 
examples include the quantum spin liquid phase of a magnetic insulator, the "strange metals" occurring in the normal 
state of the high temperature cuprate superconductors, and heavy electron systems near a quantum phase transition. 
There has been an accumulation of examples, but so far no satisfactory theoretical framework exists to describe them. 
Finding solvable examples of quantum phases with no quasiparticles should provide a valuable guide for the search 
for such a framework. 

During the last decade, developments in string theory have revealed surprising and profound connections between 
gravity and many-body systems. The so-called gauge/gravity duality (or AdS/CFT, or holographic duality), relates 
a classical gravity theory in a weakly curved (d + l)-dimensional anti-de Sitter (AdS^+i) spacetime to a strongly- 
coupled d-dimensional quantum field theory living on its boundary [1-3]. This suggests that complicated questions 
in strongly interacting many-body physics can be mapped to simple problems in classical gravity, i.e. to geometry, 
potentially giving rise to a new paradigm for approaching many-body dynamics. In particular, black holes appear to 
play a universal role in characterizing quantum phases without quasiparticles, giving important insight into dynamical 
behavior of such systems. 

Black holes are known for their simplicity and universality; assuming rotational symmetry, the geometry of a black 
hole is fully specified by its mass and certain conserved charges such as electric charge, independent of other details of 
a gravity system (such as matter content). Moreover, for a given mass, there is a maximal allowed charge. A quark- 
gluon plasma has almost equal number of "quarks" and "anti-quarks" and the temperature is the main characteristic 
energy scale. Such a state is described on the gravity side by a black hole of almost zero charge, i.e. a Schwarzschild 
black hole. In contrast, systems which are of interest in condensed matter physics, like strange metals, are finite 
density systems, with temperature much smaller than the chemical potential (a strange metal has only electrons 
but no positrons!). On the gravity side these finite density systems are described by black holes having close to the 
maximum charges. Thus, in some sense black holes provide a unified description of such seemingly distant systems 
like a QGP and a strange metal. 

Using gravity dual to explore dynamical properties of holographic QGPs has provided important insights into the 
understanding of the QGP at RHIC (see e.g. [4, 5] and references therein). The goal of these notes is to bring these 
holographic tools to bear on strongly correlated systems at finite density. In particular, we review recent progress in 
finding new universal classes of non-Fermi liquids and novel quantum phase transitions using this approach. 



One of the cornerstones of condensed matter physics is Landau's Fermi liquid theory, which underlies our under- 
standing of almost all metals, semi-conductors, superconductors and superfluids. Since 1980's there has been an 
accumulation of metallic materials whose thermodynamic and transport properties differ significantly from those pre- 
dicted by Fermi liquid theory. These so-called non- Fermi liquids, which include the strange metal phase of cuprate 
superconductors and heavy fermion materials near a quantum phase transition, present important theoretical chal- 
lenges. Here we first highlight some salient aspects of the Fermi liquid theory and then point out where it breaks 
down. 

Let us start by recalling the story for a non-interacting Fermi gas (e.g. a gas of electrons in a box) for which the 
many-particle states can be obtained by simply filing single-particle energy eigenstates following the Pauli exclusion 
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principle. The ground state is then given by filling all the (single-particle) states inside a sphere 1 in momentum space 
with radius k F determined by the density of fermions and with all states outside the sphere empty. The locus of 
points in momentum space at the boundary of this sphere, k = \k\ = k F: is called the Fermi surface. The low-energy 
excitations of the system are given by either filling a state slightly outside the Fermi surface or removing a fermion 
from a filled state slightly inside the Fermi surface, and are called a particle and hole respectively. These excitations 
are gapless and have a linear dispersion (for k — k F <C k F )\ 

e(k) = E (k) -fi = — {k - k F ) = v F (k - k F ) (1.1) 
rn 

where Eo(k) = denotes the free particle energy and we are working in the grand canonical ensemble, in which the 

energy of a excitation is measured from the chemical potential \i = E F = Particles and holes are distinguished 
by the sign of k — k F . At a field theoretical level, these excitations manifest themselves as poles in the complex 
frequency plane of the retarded Green's function Gr(uj, k) for the electron operator in momentum space 

G *M) = 1 (1.2) 

uj — e[k) + z0 + 

where uj is again measured with respect to the chemical potential \i. Note that the retarded function is relevant as it 
describes the causal response of the system if we "add" an electron to the system. Fourier transforming (1.2) back in 
time we see that the propagator describes the propagation of a particle of energy e(fc): 

G R (t,k) ~ OWe-'W* (1.3) 

which is of course as it should be for a free theory. 

The situation becomes complicated once we switch on interactions between fermions, since the notion of single- 
particle states no longer makes sense. While one may intuitively expect that the qualitative picture for a non- 
interacting gas should still apply for sufficiently weak interactions, it is a priori not clear what should happen at finite 
or strong couplings. Landau's Fermi liquid theory postulates that the above qualitative picture for a non-interacting 
gas in fact persists for generic interacting fermionic systems. In particular, it assumes that 

1. The ground state of an interacting fermionic system is characterized by a Fermi surface in momentum space at 
k = k F . 

2. Despite (possibly strong) interactions among fundamental fermions, the low energy excitations near the Fermi 
surface nevertheless behave like weakly interacting particles and holes, which are called quasi-particles and 
quasi-holes. They have the same charge as fundamental fermions and satisfy Fermi statistics. The dispersion of 
a quasi-particle (similarly for a quasi-hole) resembles (1.1) in free theory 

e(k) = v F (k - k F ) + . . . v F = — , (1.4) 

where ra* can be considered as the effective mass of the quasi-particle and is in general different from the original 
fermion mass m, from renormalization by many-body interactions. 

Given the existence of a Fermi surface, the second postulate above is self- consistent. When turning on interactions, 
a particle (or hole) can now decay into another particle plus a number of particle- hole pairs. But it is not difficult 
to check that the exclusion principle constrains the phase space around a Fermi surface so much, that given generic 
local interactions among particles and holes, the decay (or scattering) rate of a particle (or hole) obeys 

e 2 

T <Ce (1.5) 

[i 

where e is the energy of a particle (or hole). Thus, despite potentially strong interactions, particle or hole excitations 
near the Fermi surface have a long lifetime and an approximate particle picture still applies. (Below for simplicity 



1 In these lectures we will only consider rotationally invariant theories, hence we will not encounter Fermi surfaces with shapes other than 
the sphere. 

2 Note that while all excitations have positive energies, it is convenient to represent a hole using negative energy, i.e. in terms of a filled 
particle state. 
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we will refer to quasi-particles and quasi- holes collectively as quasi-particles.) Thus for a Fermi liquid, equation (1.3) 
should be modified to 

G R (t, k) ~ e-^-S* , T ~ e 2 (k) (1.6) 
which implies that near the Fermi surface the retarded function for the electron operator should have the form 

G R (u>, k) = f-TZw M + ' ' ' ' ( L7 ) 



with the self-energy E(cj, fc) 



S = £ ~ ic 2 . (1.8) 



The residue Z < 1 of the pole, which is called the quasiparticle weight, can be interpreted as the overlap between the 
(approximate) one-quasiparticle state with the state generated by acting the electron operator on the vacuum. 

The concept of quasi-particles is extremely powerful and makes it possible to develop a general low energy theory 
- Fermi liquid theory - independently of the precise microscopic details of a system. With some phenomenological 
input, the theory can then be used to predict essentially all the low energy behavior of the system. For example, the 
theory predicts that the specific heat is linear in temperature (see e.g. [6]) 

C e = 7 T + ... 7~m* (1.9) 

and that the low temperature resistivity increases with temperature quadratically 

p e =p + AT 2 + ... . (1.10) 

The theory has been tremendously successful in explaining almost all metallic states in nature. 

It is important to emphasize that Fermi liquid theory does not require interactions among the fundamental con- 
stituents to be weak. For example for 3 He, one finds that ra* = 2.8ra# e , indicating that interactions among 3 He 
atoms are clearly not weak. There also exist so-called heavy electron compounds for which the effective mass for 
electron quasi-particles can be as large as 10 2 — 10 3 times of the electron mass. 

It is rather remarkable that weakly interacting quasiparticles can emerge as the low energy collective excitations 
of a strongly interacting many-body system, with only certain parameters (such as the effective mass) renormalized 
compared to the fundamental constituents of the system. The self-consistency of Fermi liquid theory can also be 
understood from an effective field theory perspective using the renormalization group [7-9]. Assuming the existence 
of quasi-particles, one can then try to write down the most general local effective field theory for them. One finds 
that due to kinematical constraints from the Fermi surface all interactions are irrelevant at low energies except 
for the forward scatterings 3 and BCS-type pairing instabilities leading to a superconductor. Note that while the 
renormalization group analysis shows that the Fermi liquid theory is a stable fixed point (up to superconducting 
instabilities), it does not tell us whether or why a specific microscopic theory will flow to this fixed point. 




FIG. 1. Left: a cartoon picture of the phase diagram of cuprate superconductors; Right: a cartoon picture of linear temperature 
dependence of the resistivity in the strange metal phase. 



3 The forward scatterings give rise to interactions among quasi-particles via their densities, which are incorporated in the Fermi liquid 
theory. Note that such interactions do not give rise to widths. 
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Even though the Fermi liquid theory has been tremendously successful, nature always has interesting surprises 
for us. Since the 1980's, there has been an accumulation of metallic materials whose thermodynamic and transport 
properties differ significantly from those predicted by Fermi liquid theory (see e.g. [10]). One prominent class of 
examples of these so-called non- Fermi liquids is the strange metal phase of the cuprate superconductors, which refers 
to the metallic state above the superconducting transition temperature T c near optimal doping (see Fig. 1). The 
strange metal phase exhibits thermodynamic and transport behavior significantly different from those of an ordinary 
metal. A particularly striking property of the strange metal phase is that the electrical resistivity increases linearly 
with temperature, in contrast to the quadratic temperature dependence (1.10) of an ordinary metal. This remarkably 
simple behavior is very robust, existing over a wide range of temperatures (sometimes until a material melts), and 
universal, appearing in all cuprate superconductors. Nonetheless, it has resisted a satisfactory explanation for more 
than twenty years. 

The anomalous behavior of a strange metal implies that the Fermi liquid theory breaks down. The immediate 
question is whether one or both of the postulates stated earlier break down. Fortunately, this question can be 
answered with the help of photoemission experiments which can probe directly a Fermi surface and its low energy 
excitations. In an Angular Resolved Photoemission Spectropy (ARPES) experiment, incident photons knock out 
electrons from the sample and the intensity I(uj,k) of the electron beam is proportional to A(uj , k) f (uj , k) , where 
f(uj,k) is the Fermi-Dirac distribution and A(u,k) is the electron spectral function defined by 

A(u,k) = -lmG R (uj,k) . (1.11) 

7T 

Near the Fermi surface, the spectral function for a Fermi liquid (1.7) is given by a Lorentzian peak centered at 
e(k) = vp(k — kp) with a width proportional to e 2 . In particular, as the Fermi surface is approached 

A(uj,k) ZS(uj -v F (k-k F )), (k k F ) . (1.12) 

For a non- Fermi liquid, it is a priori not clear what to expect of A(lj, k). Nevertheless, as a working definition, a 
Fermi surface k = kp can be defined as the surface in momentum space where there are gapless excitations, which in 
turn should result in non-analytic behavior of A(uj,k) at k = kp and uj = 0. Thus a Fermi surface and low energy 
excitations around it can be diagnosed by examining the non-analytic behavior of A(uj, k) obtained from ARPES 
experiments. 

For high T c cuprates in the strange metal region, ARPES experiments indicate that a Fermi surface still exists, but 
excitations exhibit a much broader peak than that for a Fermi liquid. The experimental results can be fit well to the 
following expression [11], postulated as "Marginal Fermi liquid" (MFL) in [12], 

uj — vp(k — hp) + E(u;, k) 

with the self-energy E(cj, k) given by 

Ti{uj) ~ cuj \oguj + duo, c real, d complex. (1-14) 

We see from (1.13) that the system appears to possess gapless excitations of dispersion relation uj = vp(k — kp). 
However, the decay rate Y of such excitations, which is given by the imaginary part of E, is now proportional to uj in 
contrast to uj 2 (1.8) for a Fermi liquid. The decay rate, which is comparable to uj, is so large, that an excitation will 
already have decayed before it can propagate far enough (i.e. one wavelength) to show its particle-like properties. As 
a result, such an excitation can no longer be treated as a quasiparticle. Also note that the residue for the pole in the 
complex plane scales like Z ~ l/log(fc — kp) as the Fermi surface is approached and thus the quasiparticle weight 
vanishes logarithmically. Thus at the Fermi surface the overlap of an excitation with original electrons vanishes. 
Mathematically, the singularity of A(u, k) at k = kp and uj = is much softer than that (1.12) for a Fermi liquid. 

Thus the strange metal phase of cuprates has a Fermi surface (i.e. there still exist a surface in momentum space 
which has gapless excitations), but the quasiparticle picture breaks down. In other words, we are dealing with a Fermi 
surface without quasiparticles. From the perspective of renormalization group, it must be that the system is flowing 
to a nontrivial fixed point distinct from that for a Fermi liquid. 

Known field theoretical examples of non-Fermi liquids which also exhibit the phenomenon of a Fermi surface 
without long-lived quasiparticles include Luttinger liquids in 1 + 1 dimensions, and a free fermion gas coupled to some 
gapless bosonic excitations, which can be either a transverse gauge field or certain order parameter fluctuations near 
a quantum critical point (see e.g. [13-16]). In the latter case gapless bosonic excitations generate nonlocal interactions 
among fermionic excitations near a Fermi surface and destroy the quasiparticles. Neither of these, however, is able 
to explain the behavior of a strange metal phase. The former is specific to 1 + 1 dimensional kinematics. In the 
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latter class of examples, the influence of gapless bosons is mostly along the forward direction, and is not enough, for 
example, to account for the linear temperature dependence of the resistivity. Even for this class of examples, there 
still exists significant technical challenges in understanding the precise nature of the IR fixed point [14-16]. 

B. Challenge: novel quantum phase transitions 

We now switch gears, turning to another challenge faced by contemporary condensed matter theory. Consider a 
continuous phase transition; these are generally be driven by fluctuations that can be both quantum and thermal. 
Thermal phase transitions are by now very well developed, based on the so-called Landau- Ginsburg- Wilson paradigm, 
which consists of the following postulates: 

• Different orders are characterized by different symmetries. Phase transitions between orders is a consequence 
of symmetry breaking. 

• The symmetry breaking is characterized by an order parameter and critical behavior for the phase transition is 
solely captured by low energy dynamics of the order parameter. For example, the free energy can be obtained 
from 




(1.15) 



where S[<j>] is the effective action for the order parameter (j). 4 

Continuous quantum phase transitions, which happen at zero temperature from tuning non-thermal control param- 
eters, have also been traditionally formulated within the Landau paradigm [17-24]. Now the Euclidean time direction 
is uncompact and the quantum critical fluctuations of the order parameter live both in space and in time. As the 
critical point is approached in addition to a divergent correlation length £, one also finds a divergent correlation 
time £ r ~ with z the so-called the dynamical exponent. Thus the order parameter now lives effectively in d + z 
dimensions, where d is the number of spatial dimensions and z is the effective time dimension. Again one expects 
that the critical behavior be captured by a path integral like that in (1.15) with now (j) depending on both space and 
time. 

This paradigm, however, has been challenged on several fronts in recent times. Firstly, there are quantum phases 
of matter whose order cannot be captured by symmetries alone, with examples including topologically ordered states 
such as the quantum Hall phases [25, 26]. Phase transitions among or out of such phases cannot be described in 
terms of order parameter fluctuations. Secondly, during the last decade, heavy fermion compounds have provided a 
variety of fascinating examples of quantum phase transitions (see e.g. [23, 27, 28] for recent reviews), which appear 
not to be explainable using the standard Landau paradigm. In these examples, while the phase transition can still 
be characterized as symmetry breaking, new modes, which are inherently quantum and are beyond order-parameter 
fluctuations, appear to emerge as part of the quantum critical excitations. 5 These experimental examples also inspired 
new theoretical findings of quantum critical points lying outside the Landau paradigm, including "local quantum 
critical points" [29] and "deconfined quantum critical points" [30]. 

Quantum criticality is hence considerably richer and more delicate than its thermal classical counterpart. In 
particular, quantum critical behavior has been seen to be associated with some of the most interesting phenomena in 
condensed matter physics including non-Fermi liquids, novel superconductivity etc. In turn, new methods are needed 
to search for, study, and characterize strongly coupled quantum critical systems. 

For a quantum phase transition which does have an order parameter, an important set of observables to characterize 
the dynamical nature of the quantum phase transition are susceptibilities of the order parameter. Suppose the order 
parameter is given by the expectation value of some bosonic operator O, then the corresponding susceptibility x(cj, k) 
is given by the retarded function for 0, which captures the linear responses of the system to an infinitesimal source 6 
conjugate to O. 

In a stable phase in which O is uncondensed, turning on an infinitesimal source will result in an expectation value 
for O which is proportional to the source with the proportional constant given by the susceptibility. However, if the 
system has an instability to the condensation of (9, turning on an infinitesimal source will lead to modes exponentially 



4 Note that one can ignore the dependence of <j> on the time direction as such dependence only give rise to gapped excitations due to 
compactness of the Euclidean time direction at a finite temperature. 

5 In some examples, the phase transition can involve a jump in the size or shape of the Fermi surfaces. 

6 For example if O is the magnetization of the system, then the corresponding source is a magnetic field. 
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growing with time. Such growing modes are reflected in the presence of singularities of x(cj, k) in the upper complex 
cj-plane. Thus the dynamical susceptibility provides an important diagnostic of potential instabilities of a system. 
Similarly, at the onset of an instability (i.e. a critical point, both thermal and quantum), the static susceptibility 
typically diverges, reflecting that the tendency of the system to develop an expectation value of O even in the absence 
of an external source. The divergence is characterized by a critical exponent 7 (see Appendix F for a review of 
definitions of other critical exponents) 

x(fc = , w = 0)~ \g-g c \-l (1.16) 
where g is the tuning parameter (which is temperature for a thermal transition) with g c the critical point. 



C. "Local" quantum critical behavior 



A particularly interesting example which illustrates the nontrivial nature of quantum phase transitions for heavy 
fermions is provided by CeCuQ- x Au x with the doping parameter x as the tuning parameter [31]. There is a quantum 
critical point at x c = 0.1 above which the system is described by an antiferromagnetic phase. One finds near the 
quantum critical point the dynamical susceptibility behaves as [31] 

k) = — 15 — (1.17) 

(k-Q)i+x c -x + T«f(%) 

where Q is the ordering vector for the antiferromagnetic order. Equation (1.17) exhibits an interesting contrast 
between momentum and frequency dependence: 

1. The spatial part of the susceptibility exhibits ordinary mean-field behavior. 

2. The self-energy depends only on frequency and exhibits nontrivial uj/T scaling (the exponent a = 0.75). 

The above behavior has been called "local quantum critical behavior" [29, 32], and cannot be explained using the 
standard Landau paradigm. 

Interestingly, local quantum critical behavior is also present in the retarded response function for the Marginal 
Fermi Liquid (1.13). The self-energy (1.14) and its finite temperature generalization depend only on uj and exhibits 
uj/T scaling [12]. In contrast, the spatial momentum dependence is the same as that for a free fermion Fermi surface. 
Note that in [12], the fermionic response function (1.13) can be obtained if we couple the electrons to spin and charge 
fluctuations whose spectral function satisfies 

TrnWnU /xo| for M<T 

ImxM ~ I XQsigncj for T <<c M <<c ^ (1.18) 

where xo is a constant and uj c is a UV cutoff. Note that this expression is /c-independent and again exhibits local 
quantum critical behavior. 

The local quantum critical behavior described above naturally arises in theoretical models of Dynamical Mean Field 
Theory (DMFT) [33] , based on a large spatial dimension mean field approximation. In DMFT one maps a lattice many 
body problem into an impurity problem and the self-energy arises from that of the impurity. Thus by construction 
the self-energy is momentum independent and its scaling behavior can be modeled from the critical behavior of the 
impurity. Note that since an impurity typically has a nonzero entropy, in general such an approximation leads to a 
finite entropy density, indicating that it describes intermediate-energy rather than ground state physics. So it is an 
interesting question whether the observed local quantum critical behavior in heavy fermion materials and cuprates 
could be due to intermediate energy effects. 7 In the course of this review we will actually be led to very similar locally 
quantum critical behavior arising from holography. 



7 See also an argument which points to this direction at the end of Sec. II F. We have more to say on this issue in the conclusion section. 
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D. Scope of this review 

Holography is the surprising statement that certain quantum field theories can be exactly equivalent to theories of 
quantum gravity that live in one extra spatial dimension. This equivalence is obviously conceptually startling, having 
far-reaching implications for our understanding of quantum gravity. Remarkably, it is also practically useful; which 
stems from the fact that it is a strong/weak duality: generally when the field theoretical side is strongly correlated, 
the gravitational description is weakly coupled and tractable. This suggests that complicated questions in strongly 
interacting many-body physics can be mapped to simple problems in classical gravity, i.e. to geometry. 

In the most well-studied examples of the duality, gravity (and associated matter fields) propagating on a weakly 
curved anti-de Sitter spacetime in d + 1 dimensions (AdS^+i) is mapped to a strongly coupled conformally invariant 
quantum field theory (CFT^) that lives at its boundary (with d dimensions). Many examples are now known in 
various spacetime dimensions, including J\f = 4 Super- Yang-Mills theory in d = 4, and ABJM theory in d = 3 [34-36]. 
These theories essentially consist of elementary bosons and fermions interacting with non-Abelian gauge fields. At 
a heuristic level one may visualize such a theory as the continuum limit of a lattice system where the number of 
degrees of freedom at each lattice site is of order 0(N 2 ), with N the rank of the gauge group. The classical gravity 
approximation in the bulk corresponds to the strong coupling regime and the large N limit. In addition to these 
theories, there also exist vastly many asymptotically- AdS vacua of string theory, each of which is believed to give rise 
to an example of the correspondence, though an explicit description of the dual field theory is not known for most 
vacua. 

The goal of these notes is to bring holographic tools to bear on the problem of strongly-interacting field theory 
states at finite density, with an attempt to make contact with the physical problems alluded to in the above. We do 
this by studying the simplest finite density system arising from the duality, obtained by turning on a nonzero chemical 
potential for a U(l) global symmetry of a CFT^, which is described on the gravity side by a charged black hole. We 
ask the following questions of such a finite density system: 

1. Does the system have a Fermi surface? What are physical properties of low energy excitations near the Fermi 
surface? 

2. Does the system have instabilities to condensation of scalar operators? What is the quantum critical behavior 
near such a critical point? 

In general for a strongly coupled system, these questions are difficult to answer. Fortunately here they can be answered 
with the help of the gravity dual. In addition to thermodynamical quantities, the main dynamical observables that 
we will use to probe the system are retarded Green's function of gauge invariant operators, which could be a scalar 
or a spinor, 

G R (t, x) = iO(t) tr (p[0(t, £), O(0)]±) , (1.19) 

where p is the density matrix characterizing the ensemble and the notation [, ]± refers to a commutator for scalars 
and an anticommutator for spinors. The retarded Green's function carries all of the real-time information about the 
response of the system and its imaginary part is proportional to the spectral function (or spectral weight) A(uo, k) = 
^IirGr(uj, fc), which captures the density of states which can be accessed by acting O on the ensemble. As discussed 
in Sec. I A and Sec. IB, the retarded function for a fermionic operator can be used to probe the existence of a Fermi 
surface and the nature of low energy excitations around the Fermi surface, while that for a bosonic order parameter 
provides a convenient diagnostic for possible instabilities, and encodes important dynamical information for quantum 
phase transitions. For real-life systems the spectral function of various operators can be directly measured; for example 
that for the electron operator by photo-emission experiments and that for spin polarizations by neutron scatterings. 

We will see that, even in this simplest context, holographic systems offer us plenty of opportunities for studying 
exotic phenomena associated with Fermi surfaces without long-lived quasiparticles and novel quantum phase transi- 
tions, some of which share intriguing parallels with cuprates and heavy fermion systems discussed earlier. Gravity 
provides us with fascinating ways to characterize them. Our discussion will be mainly based on [37-4 7]. 8 In particular, 
we will find a universal intermediate- energy phase, called "semi- local quantum liquid" in [44], which does not have 
a quasiparticle description, and is characterized by: (1) a finite spatial correlation length; (2) an infinite correlation 
time and associated nontrivial scaling behavior in the time direction; (3) a nonzero entropy density. In a semi-local 
quantum liquid, the self-energy for an operator depends on spatial momentum 9 k only through the ratio k/fi, where 



8 The application of holographic methods to condensed matter physics was initiated in [48] and other reviews include [49-56] . 

9 We will only consider rotationally symmetric systems in this review, k refers to the magnitude of momentum. 
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ji is the chemical potential. Thus in the kinematic regime where momentum variation is much smaller than the 
chemical potential (e.g. near a Fermi surface), the momentum dependence can be neglected and one finds the "local 
quantum critical behavior" described in Sec. IC. For example, we will find Fermi surfaces which exhibit the MFL 
behavior (1.13) and quantum phase transitions whose dynamical susceptibility has the behavior (1.17) or (1.18). 

For a holographic system at a nonzero chemical potential, this unstable semi-local quantum liquid phase sets in 
at an energy scale of order of the chemical potential, and then orders into other phases at lower energies, including 
superconductors, antiferromagnetic-type states, and Fermi liquids. While the precise nature of the lower energy state 
depends on the specific dynamics of the individual system, the semi-local quantum liquid arises universally from these 
lower energy phases through deconfinement (or in condensed matter language, fractionalization). 

In this review we will be mainly interested in searching for new physical scenarios and dynamical mechanisms which 
can arise for a class of holographic systems, rather than detailed phase structure or dynamical effects for a specific 
theory. Thus we will take the so-called "bottom-up" approach, i.e. we will just consider a certain type of operator 
spectrum without referring to a specific theory. The results should apply to any theory which contains the given 
spectrum. This complements the so-called "top-down" approach which starts with specific string theory embeddings 
and works out the precise phase structure for these theories. See [57-65] for some examples of the top-down approach. 

The plan of the paper is as follows. In next section we discuss various aspects of gauge/gravity duality for a finite 
density system. In particular, we introduce the notation of an intermediate-energy semi-local quantum liquid phase 
and work out its physical properties. In Sec. Ill we turn to study of holographic (non)-Fermi liquids and in Sec. IV 
possible scalar instabilities and associated quantum phase transitions. We conclude in Sec. V with a summary of main 
points and a discussion of possible implications for real-life systems. 



II. GAUGE/GRAVITY DUALITY AT A FINITE DENSITY 



In this section we discuss various aspects of gauge /gravity duality at a finite density, in preparation for our discussion 
of non- Fermi liquids and quantum phase transitions in subsequent sections. 



A. Some aspects of AdS/CFT 

In this subsection we highlight certain aspects of gauge/gravity duality that will be needed for our discussion. 
More detailed expositions on AdS/CFT are available: for example, the TASI lectures of Polchinski [66]. Other recent 
reviews include [5, 50, 67]. The reader who is already familiar with the basic ideas of holography is encouraged to 
skip this subsection, and the reader who only wants to compute things can probably skip to the results (2.10) and 
(2.11). 

The gravitational action in AdS takes the form 
where R is the AdS radius, and Newton's constant 



2k 2 1 



R d-1 N 2 



(2.2) 



is small in the large N limit of the boundary theory. The simplest solution satisfying the field equations arising from 
(2.1) is pure anti-de Sitter space, which takes the form 

ds 2 = R 2 + l^dx" ^ = R? {dz2 _ df2 + dg2) (2 3) 



Here the x^ = (£, x) run over the coordinates of the field theory, and z is the extra "holographic" coordinate in the 
gravitational description, with the AdS boundary lying at z = 0. As the simplest and most symmetric solution on 
the gravitational side this geometry represents the vacuum of the dual field theory. Indeed, the vacuum of the CFT 
should be invariant under the conformal group in ^-dimensions, which is precisely the same as the isometry group of 
AdSd+i- We draw special attention to the following isometry: 

z -> Xz -> \x» (2.4) 



which represents scaling. 
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One of the most important aspects of the duality is that the bulk spacetime geometrizes the renormalization group 
(RG) flow of the field theory, which can be seen as follows. Due to the warp factor in front of the Minkowski 
metric in (2.3), local proper energy and length scales along field theory directions in AdS are related to those in the 
dual field theory by a z-dependent rescaling. More explicitly, consider a physical process with proper energy E\ oc at 
some value of z in AdS. The bulk energy E\ oc is measured in units of the local proper time dr = ^dt and when viewed 
in terms of boundary time t, it corresponds to an excitation of boundary theory energy E given by 

E = ^E loc . (2.5) 

We thus see that physical processes in the bulk with identical proper energies but occurring at different radial positions 
correspond to different field theory processes with energies that scale as 1/z. In other words, given that typical local 
bulk energies are given by the curvature scale, i.e. E\ oc ~ a field theory process with a characteristic energy 
E is associated with a bulk process at z ~ 1/E [1, 68, 69]. This relation implies that the z can be identified as 
renormalization group scale of the boundary theory. In particular, the high-energy (UV) limit E — > oo corresponds 
to z — >• 0, i.e. to the near-boundary region, while the low-energy (IR) limit E — >• corresponds to z — )• oo, i.e. deep 
in the interior. In a conformal theory, there exist excitations of arbitrarily low energies, which in the bulk is reflected 
in the fact that the geometry (2.3) extends all the way to z —> oo. If instead we consider a boundary theory with 
a mass gap m, the corresponding bulk geometry will then end smoothly at a finite value zq ~ 1/m. Similarly, at a 
finite temperature T, which provides an effective IR cutoff, the bulk spacetime will be cut off by an event horizon at 
a finite zo ~ 1/T. 

Another important aspect of the duality is the field/operator mapping, i.e. to each (conformally primary) operator 
one can associate a field in the bulk. For illustration, we will consider a scalar operator O in the boundary theory. 
This is dual to some massive bulk scalar field (j) with mass m, whose action can, for example, be written as 



S * = -\J dd+lx ^9 (W) 2 + m V + . . .) . (2.6) 

By studying the equations of motion arising from the action above one can show that near the boundary a solution 
to the scalar wave equation has the expansion 

<j>{z -> 0, x») ~ A(x)z A ~ + B(x)z A + (2.7) 

where the two exponents satisfy 



d d 2 
A± = -±vu v v = y - + m 2 R\ (2.8) 

Then one can show that: 



• A + is the conformal dimension of the dual operator O. 

• A(x), which is the coefficient of the more dominant term in (2.7) (in the z — » limit) can be identified as the 
source for (9, i.e. a nonzero A(x) corresponds to adding to the boundary theory Lagrangian a source term 

^boundary = J d d xA(x)0(x) . (2.9) 



• B(x), which is the coefficient of the subdominant term in (2.7) (in the z — >• limit), can be identified with the 
expectation value of the operator 0(x), i.e. 

(0(x)) = 2vuB(x). (2.10) 

For example, if one finds a regular solution where A = but 5^0, this implies that the operator O has 
spontaneously developed an expectation value, in the absence of a source. 

• Given (2.9) and (2.10), one immediately sees that the linear response function (in momentum space) for O 
should be given by 



(2.11) 



12 



where B(lj, k) and A(u, k) are the Fourier transform along the boundary directions of the corresponding quantity 
in (2.7). In order to determine the ratio in (2.11), one needs to provide an additional boundary condition at 
the interior of the spacetime, typically a regularity condition. In the situation that a horizon develops in the 
interior, the boundary condition for computing a retarded correlator is that the solution for <\> should be inf ailing 
at the horizon [70]. 

• For vjj G (0, 1), both terms in (2.7) are normalizable and there are two ways to quantize <j) by imposing Dirichlet 
or Neumann conditions at the AdS boundary, which are often called standard and alternative quantizations 
respectively. This means that depending on the boundary conditions the exact same bulk gravity action can 
be dual to two different CFTs. Our discussion above refers to the standard quantization. In the alternative 
quantization the roles of A(x) and B(x) in (2.7) are reversed, i.e. B(x) corresponds to the source and A(x) to 
the expectation value. Also the conformal dimension dimension for O is now given by A_ which lies between 
(| — 1, f ); note that the lower limit (corresponding to vjj — » 1) is now approaching that of a free massless 
scalar in d spacetime dimension. In the alternative quantization, the double trace operator O 2 is relevant, upon 
turning on which the system flows to the standard quantization in the IR. Thus this construction provides us 
with a relevant operator that can be tuned; this will play an important role in our discussion of quantum phase 
transitions below. See Appendix C for more discussion. 

At a practical level, the above (and their generalizations to spinor fields) are all that we will need to perform the 
computations in the rest of these lectures. Note that while the factors of 2vjj in (2.10) and (2.11) can be important, 
we will generally omit them in our discussion below for notational simplicity. Before moving on we summarize our 
understanding so far: 

Boundary CFT d Bulk AdS d+i 

Complicated many-body dynamics Physics encoded in geometry 

RG flow of the system Bulk "holographic" coordinate z 

Bulk scalar field <fi <^ CFT scalar operator O 

N — ^ oo and strong coupling classical GR 

The last entry in the table deserves further explanation; generally the specific field theories we are considering are 
gauge theories with gauge group SU(N). It is a well-known fact [71] that such theories are in some sense classical 
(but not trivial!) in the large- N limit, by which we mean that correlators factorize and that the physics should be 
determined by some classical equations of motion, which are not those of the naive h —> classical limit. AdS/CFT 
now tells us that the relevant classical equations of motion are actually those of higher-dimensional Einstein gravity. 
One must always keep in mind the important caveat when comparing to more-conventional physical systems that 
the observables computed using classical AdS/CFT are generally exact only in the N — > oo limit. Indeed there exist 
interesting effects that are visible only at finite TV; these are dual to quantum effects in the bulk, and we will discuss 
some of them in later sections (see e.g. Sec. HID for perturbative corrections and Sec. HIE for non-perturbative 
effects). 

B. Finite density states and AdS2 

Our primary interest is to study systems at a finite density (and temperature), the simplest way of realizing this in 
AdS/CFT is to consider a conformally invariant theory with a 17(1) global symmetry and turn on a nonzero chemical 
potential \i for the £7(1). The corresponding conserved current J fJ, (x) is dual to a bulk gauge field Am(z,x), and the 
global [7(1) symmetry in the field theory is represented by a [7(1) gauge symmetry in the gravitational description. 

Turning on a finite chemical potential \i for the current corresponds to perturbing the boundary theory by the 
operator 

^boundary = fJL J d d X J* . (2.12) 

Being conserved, J 1 has scaling dimension d — 1 and is thus a relevant operator. As a result, under (2.12), the 
system should flow nontrivially in the IR. Note that since (2.12) also breaks Lorentz symmetry, one expects the flow 
to be different in the time and spatial directions. For a strongly interacting many-body system, understanding the 
IR physics under the relevant deformation (2.12) is very difficult, and in general is not possible using conventional 
techniques. Fortunately, for a theory with a gravity dual, we can readily extract answers from the gravity side, to 
which we now turn. 
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Analogously to the scalar case discussed in the previous subsection, the boundary value of the gauge field is 
equal to the value of the field theory source for J M : thus the perturbation (2.12) means that we should study a classical 
gravity solution where the gauge field A t (z) at the boundary takes some nonzero value A t (z — » 0) = \i. We will focus 
on the minimal sector on the gravity side which is (2.1) but now with the addition of the Maxwell term for the gauge 
field, 



-9 



_ d(d-l) R 2 

7v H —3 1 n-^MN-r 

R 2 9f 



MN 



(2.13) 



^i? is a bulk gauge coupling; the factors of R in its definition combine with the dimensional Newton's constant k? out 
front to make it dimensionless. The solution that satisfies the relevant boundary condition is the Reissner-Nordstrom 
charged black hole 10 [72, 73], whose metric and background gauge field are given by 



R\ ,,.0 . ^ . tfdz 2 
z 2 f 



ds 2 = ^(-fdt 2 +dx 2 ) + (2.14) 

with 



f=l + Q 2 z 2d - 2 -Mz d , A t = fi (l--5=2 > ( 2 - 15 ) 

V Zq J 

where Q, M, zo are constants. 

In a specific string realization of such a gravity dual, there are many other matters fields in addition to those 
in (2.13) and in general it is not warranted that (2.14) will lift to a solution of the full string theory. 11 Nevertheless, 
there are many examples (including certain U(l) subgroups of M = 4 SYM and ABJM theories, and an infinite 
family of theories discussed in [87]) in which (2.14) does survive the lift to the full gravity description. Within these 
examples, the charged black hole (2.14) thus provides a "universal" geometric description of many different systems 
at a finite chemical potential, independent of specific microscopic details. Different systems will have different matter 
field spectrum on the gravity side, but the underlying geometry is the same. 

Now we proceed to describe properties of (2.14). Compare it to (2.3): the nontrivial function f(z) indicates that 
the physics is changing with scale. The horizon lies at z = zq where 

f(z ) = 0, M = Zq 4 + Q 2 Zq~ 2 . (2.16) 

Q can be expressed in terms of the chemical potential n as 



2{d - 2) n 
d-1 9Fzt 2 ' 



(2.17) 



The temperature of the system can be identified with the Hawking temperature of the black hole, which can be written 

as 

T =M i - d -^ Q ^)- (2 - i8) 

The solution (2.14)-(2.15) has two independent parameters \i and zo, with \i setting the unit of scale and the solution 
is characterized by the dimensionless number /jlzo (which can in turn be expressed in terms of the dimensionless ratio 

The field theory quantities dual to this black hole are given by the following: 

ffd—l q 

Charge density: p = J2(d-l)(d- 2) — ^ — (2.19) 

k a g F 

J^d—l 2 

Entropy density: s = 2tt — r d^i (2.20) 

« z 

^ 2 J^d—l 

Energy density: e = — t—M . (2.21) 



10 Here we say "black hole" even though the horizon is planar with topology R d_1 ; it would perhaps be more accurate to call this solution 
"black brane." 

11 This depends on whether the Maxwell-Einstein action (2.13) is a consistent truncation of the full theory. For example, if the Maxwell 
field is coupled to a scalar field <j> as e^F 2 , then (2.14) will no longer be a solution as <fi will also be excited [97, 98]. 

12 Note that since the system is scale invariant in the vacuum, at finite density the system is characterized by the dimensionless ratio — . 
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From (2.2), all these quantities are of order 0(N 2 ). It can be readily checked from the above equations that the first 
law of thermodynamics is satisfied 

de = Tds + fidp . (2.22) 

Up to some constant factors Q is essentially the charge density of the boundary theory. It is convenient to parameterize 
it as 



by introducing a length scale z* (which is fixed by boundary charge density). In terms of z*, various quantities can 
be written as 

B*- 1 1 1 d(d-l)lfz \ d - 2 „ d ( zl d ~ 2 \ , A , 

»=^r^ o'^TAi) T -*a[ 1 -i^) <2 - 24) 

where we have introduced 

e d = - 9F . (2.25) 

y/2d(d - 1) 

Note that z* > zo as otherwise the metric will have a naked singularity (rather than a black hole). 
Prom (2.24), at zero temperature, we have zq = z* and the black hole metric can be written as 

d z 2d ~ 2 2(d-l)z d „ ( z d ~ 2 \ 



z^ 



where /j, and z* are related by 



Ai , = l=2(d-2)e 4i 4- ( 2 - 27 ) 



For convenience we introduce the appropriately rescaled /i* , which will be used often below as it avoids having various 
numerical factors and gp flying around. The corresponding field theory quantities can now be written as 

_ R^ 1 1 1 __ R^ 1 2tt _ Rd- 1 (d-1) 2 1 

Note that the entropy density is nonzero even at T = 0. This finite ground state entropy density appears to violate 
the Third Law of Thermodynamics; we will discuss the implications of this in Sec. II F. 

For small z <C z* we have / ~ 1, and the metric resembles that of pure AdS (2.3); this is telling us that physics at 
energy scales much larger than the chemical potential \i is simply that of the conformal vacuum, as expected. Let us 
first consider T = 0. At z ~ z* the geometry is very different: / has a double zero at the horizon z = z*, with 

f{z)^d{d-l) {z *~ 2 Z)2 + ..., (z^z t ). (2.29) 

This implies that the horizon is actually an infinite proper distance away, as one can readily check from (2.14). In 
fact, the near- horizon geometry factorizes into AdS2 x 



_J(_^2 + dc 2 ) + /i 2 jR 2 df2) A= Zd 

C 2 1 c 



ds z = -Jr{-dt z + dC) + ntR z dx z , A = -fdt. (2.30) 



Here we have defined a new radial coordinate £ and R2 is the curvature radius of AdS2, 

c = d (d-i)U- z y = jdjd^Tj ' (2 ' 31) 

The metric (2.30) applies to the region <C 1 which translates into fi( ^> 1. From (2.27) and (2.31), a relation 
which will useful later is that 
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See Fig. 2 for a cartoon picture of the geometry of an extremal charged black hole. 

At finite temperature, in the regime z *~ z ° <c 1, i.e. T//jl <C 1, the near horizon metric is obtained from replacing 
the AdS 2 factor in (2.30) by an AdS 2 black hole [40], 

ds 2 = § (- (l - dt 2 + ) + ^R 2 dx 2 (2.33) 

with 



c 2 V CoV i-£ 

^0 



and the temperature becomes 



2^Co 

At a temperature comparable to /i, the AdS 2 structure will be lost. 



C V Co/ d(d- 1)(^ - zq) 



(2.35) 




FIG. 2. A cartoon picture of the geometry of an extremal charged black hole. The horizon lies at an infinite proper distance 
away and in the near horizon region the warp for the spatial part approaches a constant, while that for the time direction 
shrinks to zero at the horizon. 

As discussed in the last subsection, the radial direction z can be identified as the RG scale of the boundary theory. 
The emergence of the near-horizon AdS 2 region thus holds the key to understanding the low energy physics of the 
boundary theory at finite density. The warp factor in (2.30) and the fact that the horizon lies an infinite proper 
distance away imply that the system should again possess gapless degrees of freedom. In particular, the metric (2.30) 
has a scaling isometry 

t->\t, C^AC, x^rx (2.36) 

under which only the time coordinate scales (compare this to the relevant scaling symmetry for Lorentz-invariant 
AdSd + i (2.4) in which the spatial and time coordinates scale in the same way.) Note that the AdS 2 region is in fact 
invariant under a full SX(2,R) symmetry group, and likely a Virasoro symmetry. 
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Thus gravity tells us that at low energies the boundary system flows to a fixed point which is dual to AdS 2 xR d_1 [40] . 
See Fig. 3 and 4 for illustration. This fixed point has previously been called the "IR CFTi" to indicate that there 
is conformal symmetry only in one dimension (the time direction). However, this hides the fact that the dependence 
on the spatial directions can have important consequences, and thus we believe that a more descriptive name for this 
fixed point is a semi-local quantum liquid (or SLQL for short) [44], for reasons that we will elaborate on in the next 
subsection. We will use this term in what follows. 

In addition to translational and rotational symmetries of the boundary spatial directions, the SLQL is thus char- 
acterized by conformal symmetries of a (0 + l)-dimensional conformal quantum mechanics including the scaling 
symmetry in the time direction. It is crucial to note that this conformal invariance has nothing to do with the CFT^ 
in the ultraviolet, whose conformal invariance was completely broken by the chemical potential. This new conformal 
symmetry is emergent, and has to do with the collective motion of the large number of charged excitations sustaining 
the background charge density. 

Note that the SLQL is specified by a single scale /i*. While /i* corresponds to the chemical potential for the full 
UV theory, near the IR fixed point SLQL, ji* simply defines a unit for spatial momenta. 



IR UV CFTd 




FIG. 3. Charged black hole with nonzero gauge field is dual to field theory state with nonzero chemical potential; in the infrared 
there is an emergent conformal symmetry corresponding to the AdS2 x R d_1 part of the geometry. 



AdS d+1 < ► CFT d 



AdS 2 x R**" 1 * ► SLQL 



FIG. 4. At a finite chemical potential, a CFT^ flows in the IR to SLQL. On the gravity side this is realized geometrically via 
the flow of the AdS^+i near the boundary to AdS2 x R d_1 near the horizon. 



C. Semi-local quantum liquids 

Consider an operator 0(t,x) in the boundary theory, which is dual to some bulk field <fr(t,x,z) with mass m and 
charge q. At the UV, i.e. CFT^, the dimension of O is given in terms of bulk quantities by (2.8). At the IR fixed 
point, O should match to some scaling operator $ (under scaling symmetry (2.36)) in the SLQL dual to the bulk field 
(j) in AdS 2 x R d_1 . 
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Now in AdS2 x ^ d 1 , the Fourier transform 13 fc, z) of x, z) along the spatial directions has an effective AdS2 
mass depending on spatial momentum k = |fe| (as in a Kaluza-Klein reduction), which will result in a /c-dependent 
scaling dimension for boundary operators. It is then convenient to label operators by their spatial momentum, with 
fc, z) in AdS2 x dual to some operator in the SLQL. The Fourier transform 0%(t) of 0(t, x) should then 

match to at the IR fixed point. Note that the distinction between Or and is important; (9 is an operator in 

the full theory, while only exists at the IR fixed point. 

The scaling dimension 5^ of and its correlation functions in SLQL can be obtained from solving the relevant 

wave equations for <j> in the AdS2 x geometry (2.30) [40, 44]. We present the results here, leaving their derivation 
to Appendix A. One finds that 

Sk=2+"k ( 2 - 37 ) 

with 



On0 9 1 k 2 R 2 1 / OTV> g 2 F q 2 did -I) k 2 /o ooN 



and 



k = \k\, q* = qe d , (2.39) 

where in the second equality of (2.38) we have used (2.32). The retarded function of in the SLQL 14 can be written 

as 

g k (uj) = c(u k )(-iuj) 2 ^ (2.40) 

where 

Let us now examine what (2.38)-(2.40) tells us regarding the physical properties of a SLQL. Firstly, the dimen- 
sion (2.38) depends on the momentum momentum k: as a result operators with larger k become less important in the 
IR. Also Vk decreases with q, i.e. an operator with larger q will have more significant IR fluctuations (given the same 
mass m). From the bulk point of view, this ^-dependence in operator dimension is a consequence of the fact that in 
AdS2, the electric field following from (2.30) blows up as we approach the AdS2 boundary. 15 

Secondly, the spectral weight, which is defined by the imaginary part of the retarded function (2.40), scales with uj 
as a power for any momentum /c, which indicates the presence of low energy excitations for all momenta 16 (although 
with a larger scaling dimension the weight will be more suppressed at larger momenta). 

Thirdly, in coordinate space, the system has an infinite correlation time, but a finite correlation length in the spatial 
directions. This is intuitively clear from the presence of (and lack of) scaling symmetry in the time (spatial) directions. 
It can be seen more explicitly by examining the Euclidean correlation function Ge(t = it,x) in coordinate space by 
Fourier transforming (2.40). For this purpose, we note that can be rewritten as 



where 



£ = — . 1 = - 1 . (2.43) 

^* Jm 2 R 2 - ^kf- + d(d ~ 1} V d ( d l )vk=W* 

Vk has a branch point at k = i£ -1 , which under the Fourier transform leads to exponential decay in spatial directions 
at large distances with a correlation length given by £ (see Appendix D). Explicit calculation yields that Ge(t,x) has 
two distinct regimes: 



To simplify notations, we will always use the same symbol to denote the Fourier transform of a field, distinguishing them only by their 
arguments. 

We stress that this is not the full two-point function of the operator O in the UV theory, but rather just the correlation function at the 

IR fixed point; the relation between these two is given in (2.56) and is discussed in Sec. II D. 

In other words, the electric field is non-renormalizable and specifies the corresponding boundary theory. 

Exercise for students: compare this result to more "ordinary" gapless systems, like a free massless boson or a Lorentz-invariant CFT 
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1. For x = \x\ < £, 

G ^x)^-^-. (2.44) 

2. For x ^> £, the correlation function decays at least exponentially as 

G e (t,x) ~ e~f . (2.45) 

Fig. 5 provides a heuristic visualization of the behavior: the system roughly separates into independent clusters of 
size of order £, with the dynamics of each cluster controlled by a conformal quantum mechanics, with a power law 
correlation (i.e. infinite relaxation time) (2.44). Domains separated by distances larger than £ are uncorrelated with 
one another. Given that the system has a nonzero entropy density, each cluster has a nonzero entropy that counts 
the number of degrees of freedom inside the cluster. 



FIG. 5. A cartoon picture: the system separates into domains of size £ ~ — . Within each domain a conformal quantum 



mechanics governs dynamics in the time direction with a power law correlation (i.e. infinite relaxation time). 

While (2.44) and (2.45) can be found by doing Fourier transforms explicitly, they can also be seen geometrically using 
a geodesic approximation to calculate field-theoretical correlation functions using the Euclidean analytic continuation 
of (2.30), 

R 2 

ds 2 = -JtUt 2 + d( 2 ) + R 2 fil dx 2 . (2.46) 

C 2 

Consider a geodesic that starts at ( —> 0, moves into AdS2 at finite £, eventually turns around, and returns to the 
boundary but at a spatial separation x and a Euclidean temporal separation of r. In the geodesic approximation the 
boundary theory Green's function in coordinate space Ge(%, t) ~ e -mL(r,x) w \ iere L(r^ x) is the proper distance along 
the geodesic. Since the metric (2.46) is just a direct product, we can simply find the distance moved in each factor 
and add them using Pythagoras. The distance moved in the R d_1 factor is R[i*x. A standard calculation shows that 
the distance moved in the AdS2 factor is 2R2 log (^) where e is an IR cutoff of the AdS2- 17 Note that since the AdS2 
is embedded in the full spacetime (2.14), the cutoff e should be provided by the scale ^. From IR/UV connection, 
an IR cutoff e of AdS2 translates into a short-distance cutoff in the SLQL, thus we will restrict to r > e. Combining 
these results we find 



G e (t, x) - exp ( -mJm 2 log 2 Q + ^R 2 x 2 



exp 




452 log2 (r) + |_) {2A7) 



with S = mi?2 and £ 1 = mRfi*. The geodesic approximation applies to mR ^> 1 and in this regime, (2.44) and (2.45) 
are indeed recovered from (2.47). 



The easiest way to understand this result is to note that its exponential must reproduce the conformal result, r 
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The fact that at low energies different points on the 



pd-l 



can be thought of as being in different disconnected 



domains with size £ ~ ^ can also be seen geometrically as follows. Consider two spacetime points on a hypersurface 
of given ( — >> oo. To see whether observers at those locations can communicate with each other we look at time-like 
geodesies in (2.30) which connect the two points. Simple calculations show that there is a maximal separation in 
R d_1 directions for two points to communicate with each other, given by 

1 7T 



Ax n 



= 7rR<2 



V*R ii*y/d(d- 1) 



(2.48) 



The first factor 7ri?2 in the first equality above is the time for a timelike geodesic (see Fig. 6) to approach the boundary 
and come back and the second factor is the effective velocity in R d_1 (see (2.30)). Equation (2.48) is consistent 
with (2.42) up to a prefactor. 



X, T 




X 



Horizon 



Horizon 



FIG. 6. Geodesies in the near horizon geometry. In both cases the horizon is the dotted line at £ — >• oo, and the boundary is 
to the right towards f — » 0. The left plot gives the geodesic for calculating the Euclidean correlation function (2.47); the right 
plot gives a timelike geodesic connecting two points on a constant-^ hypersurface. The maximal boundary time it takes for a 
timelike geodesic to approach the boundary and come back is ttR2- 

To summarize, equation (2.40) describes a disordered state in which the space factorizes into independent domains 
of correlation length £. Within each domain one has scale invariance along the time direction. This behavior is 
reminiscent of various theoretical models based on a large spatial dimension mean field approximation [33], such as 
the gapless quantum liquids of [74-76] and the "local quantum critical point" of [29] 18 . We should emphasize that 
the scaling behavior within each cluster here describes, however, not the behavior of a single site, but the collective 
behavior of a large number of sites (if one considers our systems as a continuum limit of a lattice) over a size of order £. 
It is clear that while there is nontrivial scaling only in the time direction, the scaling dimension (2.38) and correlation 
functions (2.40) depend nontrivially on k. As discussed earlier around (2.42) it is precisely this dependence that gives 
the spatial correlation length of the system. Also, while at a generic point in parameter space the dependence of v k 
and Qk on k is analytic and only through fc//i, as we will see in Sec. IV, near certain quantum critical points, the 
dependence becomes nonanalytic at k = and is important for understanding the behavior around the critical point. 
It is also important to emphasize that, despite the scaling behavior in (2.40) we are describing a phase, not a critical 
point. We thus call it a semi-local quantum liquid, or a SLQL for short. 

We now consider some generalizations of (2.38) and (2.40). At a finite temperature the background geometry is 
given by an AdS2 black hole (2.33) and equation (2.40) generalizes to [40, 42] 

5f( W )=J*"s t (^^) (2.49) 

where is a scaling function given by 

9b\yk,y = 4tt) k — I : — r=yi . . \ • 2.50 

r(2i/fc) r(± -v k - iq*) r{±-v k -iy + iq*) 



Some suggestions for a connection have been made in [77-80]. 
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For a spinor field one finds that, by studying the behavior of the Dirac equation in the AdS2 region, the dimension 
5 k of a spinor operator in SLQL is given by 



Note this is slightly different from the corresponding result for the charged scalar (2.38). Similarly, the retarded 
Green's function is given by 19 

g k (uj,k)=c(iy k )(-iuj) 2uk (2.52) 

with c(v k ) now given by 

[ h) " r(2, fe) r ( i - v k - * g .) ( m _ R2 _ ^ + Vk • (2 ' 5d) 

Similarly, at finite temperature, the spinor AdS2 Green's function equation (2.52) generalizes to 

g(T) {uj )=T^g f (v k ,^ f ) (2.54) 

with gf given by 

„ = (**r> t - ?)* - k - - Vk T ^ r(i + ^ - r ( j + >-t - % + %») (M5) 



( m _ ^ ^ 2 _ ^ + Vk T(2^) T(l - - iq*) T(\-v k -iy + iq* 



D. Low energy behavior of retarded Green functions 



Having understood the behavior of the retarded function at the IR fixed point SLQL, we now turn to the low 
energy behavior of the retarded function for O in the full theory, which we will denote as Gr(u, k) to distinguish it 
from Gk(w) at the IR fixed point. We will consider a scalar operator at T = and mention the generalization to a 
finite temperature and a spinor at the end. We are interested in the low frequency regime uj <C \i. Note that since we 
are working in the grand canonical ensemble the frequency uj in our discussion of this paper always refers to energy 
measured from the chemical potential. 20 

GR(uj,k) can be computed explicitly by solving the equation of motion for the corresponding bulk scalar field in 
the full charged black hole geometry (2.14) and following the procedure outlined in Sec. II A. Basically our task is to 
understand how the expansion coefficients A, B in (2.11) behave at low frequencies. We present only the final result 
here, 



ft = n 2uu M^, fe) + M", k)g k (u)^ 

a + (cj, k) + a_(a;, k)G k {oj)n*' 



g r ( U , k) = £*> T';; "rrr:\ + • ■ • > ( 2 - 56 ) 



relegating the details to Appendix B. The various objects appearing in this formula deserve explanation. G k (uo) is the 
retarded Green's function in the SLQL discussed in Sec. II C and is proportional to uj 2vk with v k defined in (2.38). 
a± and b± are quantities that arise from the solving the equation of motion in the region of the black hole geometry 
not too close to the horizon (called the UV region in Appendix B); the key fact here is that they are analytic in uj 
and have a smooth uj — ^ limit. See Appendix B for further discussion of their properties. Since a±,b± in (2.56) are 
analytic in cj, we can expand them as follows 

a+(cj, k) = a { °\k) + a^\k)uj + (k)uj 2 + ... (2.57) 



The following expressions are for half of the components. The expressions for other components are obtained by taking k — > —k. 
See [40, 42] for details. 

More explicitly, for an operator of charge q, the effective chemical potential is then fiq and uj is the energy measured from ^q. See also 
discussion below equation (Bl) in Appendix B. 
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and keep only the leading order term after which (2.56) becomes 

GrIuj, k) = a^ u 7ry\ 77^ \- ... . (2.58) 

Equation (2.56) contains two sets of data: "universal" data from Gk(u) that depends only on the IR fixed point (or 
the near-horizon AdS2 x region of the black hole geometry), and non-universal data a±,b± which depends on the 
rest of the black hole geometry and can be thought of as encoding the the effects of the UV degrees of freedom on the 
low energy physics. This is consistent with general expectations from a Wilsonian understanding of renormalization 
group flow, and we now show that the form of (2.56) follows simply from general considerations of a suitable low 
energy effective action. 

We expect that the low energy effective action of the system can be written as 

&eff = Sslql + Suv (2.59) 

where Sslql is the action for the IR fixed point, for which we do not have an explicit Lagrangian description, but as 
discussed in Sec. II C whose observables can be computable from gravity in AdS2 x Suv arises from integrating 

out higher energy degrees of freedom, and can be expanded in terms of scaling operators in Sslql- The part relevant 
for O can be written as 21 

Suv = lj A(fc,u,)jf + J V (k, u)*nJ_ U ~\J <i{k,oj)*l + ■■■ (2-60) 

where we have written the action in momentum space since the dimension of SLQL operator is momentum- 

dependent, and the integral signs should be understood as J = J duo J dk. We have introduced a source for 
and . . . denotes higher powers of and J. Note again the distinction between O and <£; $ is the operator at the 
IR fixed point to which O matches. Since we are only interested in two-point functions it is enough to keep Suv to 
quadratic order in <£> and J. For simplicity we have assumed to the quadratic order there is no mixing of O with other 
operators. 22 The discussion can be easily generalized to the situation with mixing. The functions A(/c,co>), rj(k,uj) and 
£(k : oj) parametrize the "UV data" coming from integrating out higher energy modes. In particular they should be 
real and have an analytic expansion in uj with a nonzero uj —> limit, i.e. when written in terms of t, the effective 
action (2.60) should have a well defined derivative expansion. 

The (Euclidean) two-point function for can be obtained by differentiating the Euclidean partition function of 
S e ff with respect to Jr. It contains two parts. The first part is simply A(cj, fc) from the first term in (2.60), which 
comes purely from the UV and is independent of any IR fields. The second part comes from the correlation function 
4>£ in the SLQL (up to a factor of rj). To the order written explicitly in (2.60), the full action (2.59) describes the 
SLQL deformed by a double trace operator with coupling £. The two-point function of in this deformed theory can 
be obtained from the general formula (C5) in Appendix C. Combining the two parts and analytically continuing to 
Lorentizian signature, we find that the retarded function Gr for the full theory can be written as 

G R (u,k) = A(k, U )+ fl k '"} - ■ (2.61) 

Equation (2.61) has exactly the same form as the gravity result (2.56) and in fact comparing them we can identify 

A = n1 vu — , i = H-^— , V 2 =™H? u - 2vk , W = a + b--a.b + . (2.62) 

(2f_|_ (2r_|_ 

The low energy effective action (2.59), (2.60) with identifications (2.62) can also be derived directly following the 
formalism of holographic Wilsonian RG introduced recently in [81, 82]. Note that while our main goal here is to 
derive the retarded function (2.56), the effective action (2.59)-(2.60) clearly has much wider applications and we will 
see some of them in later sections. 

Since here the spatial directions do not scale, the last term in equation (2.60) is an irrelevant perturbation from 
SLQL since has dimension \ J rV\ z > \. This is what we should expect; SLQL is the IR fixed point and the last 



Note the expressions below are written somewhat heuristically to exhibit the main structure. For example, J| should be interpreted as 
JrJ r and similarly with 

K —K ' k 

This corresponds to that the equation for the dual bulk field <p does not mix with other fields at linear order. 
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term in (2.60) is simply the lowest irrelevant operator taking it to the UV. Indeed from (2.58) (or equivalently (2.61)) 
we find that as uo — > 0, 

TmG R (u>, k) = Vl" u — ^J lmg k (uj) + ... ~ uo 2vk + ... (2.63) 

v u (a^ } (k)) 2 

where we have used (B19) and that a±,b± are real. As expected, up to a non-universal prefactor (which can be 
understood as a wave- function renormalization between the UV and IR operators), at low energies the spectral 
function is given by that of the IR fixed point. In particular, as already emphasized in Sec. II C (one paragraph before 
that of (2.42)), the system has gapless degrees of freedom at any spatial momentum. 

It turns out there are circumstances in which (2.63) does not apply, and the physics is in fact much richer. This 
will be the subject of the next subsection. 

We conclude this subsection with some generalizations: 

1. At a finite temperature T <C /i, equation (2.56) can be generalized by replacing Qk(w) with its finite-temperature 
generalization Q^\uo) of (2.49). In addition, there will also be analytic finite-T corrections to a±,b±, which 
can be expanded perturbatively in T. 

2. An analogous calculation can be done for a Dirac spinor field in the bulk. Some of the details are different but 
at the end of the day the formula (2.56) still applies to each of the eigenvalues of the Green's function, which is 
now a matrix in spinor space. We discuss this briefly in Sec. Ill A; see [40, 83] for details. 

E. Possible low energy physics: Fermi surfaces and scalar instabilities 

As discussed in Sec. I A and Sec. IB, the retarded function for a fermionic operator can be used to probe the 
existence of a Fermi surface and the nature of low energy excitations around the Fermi surface, while that for a 
bosonic order parameter provides a convenient diagnostic for possible instabilities, and encodes important dynamical 
information for quantum phase transitions. If (2.63) were the full story, then that would be somewhat boring, as 
there will be no Fermi surface or instability (thus no quantum phase transition). However, equation (2.63) breaks 
down when one of the following occurs and interesting new physics emerges, 

1. a+\k) could vanish at some momentum; 

2. Vk which characterizes the scaling dimension of an operator in the SLQL becomes imaginary. Note in this case 
a±,b± become complex. See (B21) and (B22) in Appendix B. 

From the RG point of view, item 1 implies that the effective action (2.60) breaks down as (2.62) implies that the 
coefficients of (2.60) become singular in the uj — >• limit at the momentum where a+\k) vanishes. When is pure 
imaginary (say equal to — iA/~), the dimension of in the SLQL is now complex | — iA^; a complex conformal 
dimension is somewhat peculiar and seems likely to indicate that we are no longer using the correct effective theory. 
Thus in either situation we can no longer conclude that the theory is described by SLQL alone in the IR. 

In the rest of this subsection we discuss possible physical interpretations of the above two possibilities assuming 
they do occur, leaving more in-depth study (including the parameter range over which they could occur) to Sec. Ill 
and IV. 



1. Inhomogeneous black hole hair: Fermi surface v.s. scalar instability 

We now examine the first possibility, i.e. a+\k) could vanish at some momentum. Recall from equation (B6) 
in Appendix B that is the leading expansion coefficient near the AdS^+i boundary of the UV region solution 

rj^ (B5). By construction rj^ is normalizable at the horizon. When a^} vanishes rj^ becomes normalizable also at 
the AdSd + i boundary and thus is now a genuine normalizable mode in the black hole geometry. One might call it 
black hole hair. Clearly, the condition that a+\k) = is not generic; it requires fine-tuning in one of the parameters 
entering the differential equation and so can only have solutions at most at some discrete values of momentum k. 
When this happens at a nonzero momentum, one finds inhomogeneous hair. The discussion also applies identically to 
the spinor case. 
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which can be further written as 



where 



and 



Now suppose that indeed vanishes at some value of k — kp ^ 0: i.e. a+\kp) = 0. Let us consider what 
happens to (2.56) near kp. First at uo = 0, we find simply a pole 

G R {. = 0, fc) « ft* - ^ \g - - + (2.64) 
d k a\'(k F ){k - k F ) 

For a spinor operator, this is a zero-frequency singularity at a shell in fc-space, and so is precisely the signature of 

a Fermi surface! Now further turning on a small frequency to near k ~ &p and keeping the first few terms in the 
relevant expansion of (2.56), we find 

b (0) (k P ) 

G R (uj,k)^ nl"" + v ' - =2— + ... (2.65) 

d k af(k F ){k - k F ) + ua^\k F ) + a ( °\k F )g kF (oj)^ ^ 

G R (co, k) « 7 7 t 77TT + • • • (2-66) 

k — kp — u; — ZjIcj) 

v F -f , hi —^77T, ( 2 -67) 

a { +\(k F ) d k af(k F ) 

= h^^- = hc(v kF ) I — \ , h = - . (2.68) 

/i* V^*/ d k ay(kp) 

Equation (2.66), which has a similar (but more general) form as the Fermi liquid (1.7), is indeed consistent with 
the behavior of low energy response near a Fermi surface. £(tj) can be interpreted as the self-energy for fermionic 
excitations near the Fermi surface and is controlled by the retarded function Q kp of the IR fixed point (i.e SLQL) 
evaluated at k = kp. Recall now that for a Landau Fermi liquid (1.8) we have £fl(^) ~ following from basic 
kinematics of interactions about the Fermi surface. Here, however, we have 

E - u 2uk F (2.69) 

with a nontrivial scaling exponent v kp . Thus we find a non-Fermi liquid; a system with a Fermi surface but whose 
excitations appear to not be Landau quasiparticles. The dispersion and decay width of small excitations near the 
Fermi surface can be extracted from the poles of (2.66). The nature of such excitations turns out to be rather different 
depending on the value of u kF: which we will discuss in detail in Sec. Ill, along with when such a kp could exist. 

The above discussion works well for a spinor operator, but how do we interpret (2.64) and (2.66) for a scalar 
operator? This turns out to be related to an important question for the spinor case: does the pole of (2.66) always 
lie in the lower complex cj-plane as is required from analyticity? Note that at k = kp, the pole lies at the origin of 
the complex cj-plane. If we move slightly away from k = kp the pole will move off the origin and into the complex 
uj plane. What will the motion of the pole be? Based on the following property of the scalar and spinor retarded 
functions (uj > 0) 

scalars : ImG R (-u) < 0, ImG R (u) > (2.70) 
spinors : lmG R (-u) > 0, lmG R (u) > , (2.71) 

one can show that [40] 

1. For a spinor the pole always lies in the lower-half complex cj-plane. 

2. For a scalar, the pole crosses from the upper half plane to the lower half plane as k is increased from below to 
above kp. Thus for a scalar there is a pole for k < kp in the upper half plane. 

We note that (2.70) and (2.71) simply follow from the definition (1.19) that the retarded function for a fermionic 
(bosonic) operator is defined to be an anticommutator (commutator). Although in the bulk we are dealing with 
classical equations of scalars and spinors and have not imposed any statistics, the self-consistency of AdS/CFT 
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implies that the classical equations for bulk scalars and spinors should encode the quantum statistics of the boundary 
theory. 

A pole for k < kp in the upper half cj-plane for a scalar gives rise to exponentially growing modes (recall the 
discussion of Sec. IB). Thus instead of finding a Fermi surface, black hole scalar hair at some kp ^ indicates an 
instability to condensing the scalar operator. In Sec. IV we will see that it is possible to vary parameters of the system 
so that kp can be smoothly tuned to zero, after which the instability disappears. The value of the external parameter 
at which kp = thus corresponds to a quantum critical point, whose critical behavior will be the focus of Sec. IV. 



2. Low energy effective theory 

Now let us turn to the low effective theory (2.59)-(2.62) to examine what happens when vanishes. From 

equation (2.62) and (2.57), when = 0, all the coefficients in the Sjjv diverge as uj — » 0, i.e. the effective action 
becomes non-local. This means we must have integrated out some gapless modes, which of course should be the small 
excitations around kp which we have just identified from the poles of (2.66). In order to have a well-defined local 
effective action, we need to put such gapless modes back to the low energy effective action. We will introduce a new 
field \I/ in the low energy theory and "un-integrate" (2.60). Clearly there is no unique way of doing this 23 and the 
simplest choice is 

Suv = -\ J K(w,fc)$! + j A(fc )W )%%-i j *_ jg Jf( w ,A)*E + j *J (2.72) 

where 



^ = A- 1 = f±, X = VK=^, « = C+ £ = ^. (2.73) 
6 + o + K b+ 

Now all coefficients in (2.72) have local expansion in uj. This discussion applies to both a spinor or a scalar operator. 
Near fcp, K has the expansion 

K = (uj - v F (k - k F ) + . . .) (2.74) 

All 

and thus for a spinor operator, the third term in (2.72) can simply be interpreted as the action for a free fermion near 
its Fermi surface at kp. The full low energy effective action (2.59) can now be written as [40, 43] 

Seff=S S LQL[$\+ j A(fc,w)%*g + 5f re efermion[*] + j (2-75) 

m = J dtdkV_^ (id t - v F (k - k F ) + . . .) ttg (2.76) 



where 



*->free fermion 

is the action for a free fermion near a Fermi surface and 

SSLQL = SsLQL ~\ J k(w, (2.77) 

is the action for $ which is SLQL deformed by a double-trace term. The low energy action (2.75) describes a free 
fermion field ^ hybridized with a strongly coupled fermionic operator $ whose dynamics is in turn controlled by SLQL. 
(2.75) also offers a simple field theoretical interpretation for the expression (2.66): the self-energy simply arises 
from the mixing of ^ with $ and is thus given by the retarded function of see Fig. 7. As we will see in Sec. Ill, the 
action (2.75) also provides an efficient way to understand the behavior of small excitations near the Fermi surface. 

The above discussion also applies to the scalar case. For kp ^ 0, one is, however, dealing with an unstable state 
and its physical meaning is not so clear. But for kp = 0, i.e. near a quantum critical point mentioned earlier, then 
^ can be interpreted as the order parameter and the analogue of (2.75) provides an effective action for the order 
parameter. In particular, the free-fermion part in (2.75) becomes the standard Landau-Ginsburg action. We will 
discuss this further in Sec. IV. 



We can for example make a field redefinition in tp as ip — > + 
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Gr = 



~r — (slql) |- 




FIG. 7. The action (2.75) gives a simple field theoretical interpretation of (2.66) through hybridization between a free and a 
strongly interacting sectors. The propagator for the free fermion is given by K (2.74). 

3. Pair production by a black hole: complex IR dimensions 

Let us now consider what happens when the conformal dimension \ + of an operator in an SLQL becomes 
complex. Recall that for a scalar the conformal dimension is given by (2.38), which we copy here for convenience, 



gjg 2 , d(d-l) | k 2 



2 R 2 



with a similar expression (2.51) for a spinor. We thus find that for 



> 



Vk becomes imaginary for k 2 < k 2 with 



Vk = — i\k, A/c = 



-VW^k 2 . 



(2.78) 



(2.79) 



(2.80) 



y/d(d- 

Note on the gravity side, equation (2.79) corresponds to violating the Breitenlohner-Freedman bound [84] of AdS2 [85- 
88]. 

For a spinor (2.79) is replaced by 



fc 2^(4^-m 2 i? 2 ') >0. 



(2.81) 



Basically we see that the background electric field acts through the charge as an effective "negative mass square", 
making it possible to have a negative total mass square and resulting in an imaginary conformal dimension. 
Using (2.80) we now find from (2.58) the leading small frequency behavior 



G R {LO,k)=^ 



-2i\k 



-2iA fc * 



(2.82) 



As discussed in Appendix B, for imaginary v^, a± \ b±^ also become complex and in fact 

af = (a^Y bf = (b™r ■ 
As a result, even at uo = we find that the spectral weight is non- vanishing for k < k Q 

& (0) 

Im G r (uj = 0, k) = tfT u lm -±y + . 



(2.83) 



(2.84) 



We stress that this result - a finite spectral density at precisely zero frequency - is somewhat unusual, and its 
interpretation will be discussed further in Section HIE. At finite frequency (2.82) is oscillatory; it is periodic in logo; 
with a period given by = ^. In particular, it is invariant under a discrete scale transformation 



(2.85) 
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We will refer to the region k < k Q as the oscillatory region. 

It can be readily checked that (2.82) has an infinite number of exponentially separated poles in the complex cj-half- 
plane; as (2.85) shows, there is an accumulation point at uj = 0. Again using (2.70) and (2.71) one finds that for a 
spinor all the poles lie in the lower half cj-plane while for a scalar they all lie in the upper half cj-plane [40] . Thus we 
immediately conclude that for a scalar, in the parameter region (2.79), there is an instability to the condensation of 
the scalar operator. While there appears no obvious pathology for a spinor operator at the current level of discussion, 
as we will see below and in Sec. Ill, there is nevertheless an instability, and at an exponentially small (in terms of TV) 
scale, the system becomes a Fermi liquid. 

It is also instructive to understand the above phenomena for both fermions and bosons in terms of bulk language. 
For this purpose it is convenient to slightly rewrite the AdS2 x metric (2.30) by introducing 



C = R 2 e y (2.86) 



after which (2.30) becomes 



ds 2 = -e 2y dt 2 + R\dy 2 + [i\R 2 dx 2 , A = ^e~ y dt, . (2.87) 

We note that the spatial part of (2.87) is flat and the chemical potential hl for a local observer with charge q in the 
bulk is given by 

\Q\e d g F \q\ {2M) 



R2 V2R 

where in the last equality we have used (2.32). Note that [1l does not depend on the radial direction y of the AdS2- 
For a fermion, in the Thomas-Fermi approximation, we fill all the states up to hl which is the Fermi energy. This 
gives the local Fermi momentum 



where the local bulk momentum k^ is related to k of the boundary and radial momentum k y conjugate to y as 

*>=^ + I- <2 - 9o) 

Note that the bulk Fermi surface defined by (2.89) is one dimensional higher than a momentum shell in the boundary 
theory. Thus projecting it to the boundary we conclude there should be gapless excitations for all k satisfying 



/ 2 2 

k<fi*J^--m 2 R 2 . (2.91) 

See Fig. 8. Note that the right hand side of the above equation is precisely k Q introduced in (2.81). This explains the 
nonvanishing of spectral weight in (2.84) as a consequence of the projection of a bulk Fermi surface to the boundary 
theory. The oscillatory behavior in (2.82) can similarly be understood in terms of the bulk Fermi surface when taking 
into account of that the local proper energy is red-shifted by a factor e~ y compared with the boundary frequency uj 
(see (2.87)) and the behavior of local wave function. 

For a scalar, when the local chemical potential (2.88) exceeds the mass of a charged particle 

> m (2.92) 
V2R J 

the system will simply Bose condense. Equation (2.92) agrees with condition (2.79) up to an additive constant 24 and 
gives the bulk origin of the scalar instabilities we found earlier from the poles of boundary retarded functions. 

For a charged operator, the above stories can also be understood more heuristically from the black hole geometry 
as follows. When (2.79) and (2.81) are satisfied one can show that the corresponding quanta for the bulk scalar 
(or spinor) field can be pair produced by the electric field of the black hole [91]. Suppose that the black hole is 



24 Note that for our local approximation to work well, we need the typical local proper wave length to be much smaller than the curvature 
radius R, which translates into gpq^mR ^> 1. 
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FIG. 8. In the AdS 2 x R 2 (take d = 3 for example) region of the extremal black hole geometry, at each point in the bulk there is 
a local three-dimensional Fermi surface with Fermi momentum k , which upon projection to the boundary theory would result 
in a Fermi disc, in which there are gapless excitations at each point in the interior of a disc in the two-dimensional momentum 
space. 

positively charged, then the negative-charged particle in a pair will fall into the black hole, while the positive-charged 
one moves to the boundary of the spacetime. However it cannot escape, as the curvature of AdS pulls all matter 
towards its center, and thus the particle will eventually fall back towards the black hole. It then has some probability 
of falling into the black hole or being scattered back toward the boundary. For a scalar it turns out the so-called 
superradiace is happening, i.e. for a wave incident on the black hole, the reflection amplitude becomes greater than 
the incoming amplitude. Thus this pair production process will be magnified and leads to an instability. In contrast, 
for a spinor, there is no superradiance, and this process will eventually reach an equilibrium with a positively charged 
gas of fermionic quanta hovering outside the horizon. The ground state of this fermionic gas is then described by a 
bulk Fermi surface described earlier. 

For a charged fermionic field, the induced fermionic density scales as 0(N°) in the large N limit as the pair 
production is a one- loop process in the bulk. Since the charge density carried by the black hole, given in (2.28), scales 
as 0(A^ 2 ), naively one might conclude that the fermionic backreaction can be ignored at leading order. However, as 
we will discuss in more detail in Sec. HIE, the induced local bulk density is independent of radial coordinate the 
integration of which over y then gives an infinite answer, signaling an instability [95]. In Sec. HIE we show that after 
taking into account of the fermionic backreaction, SLQL becomes a Fermi liquid at very low energies. 

4- Summary 

To summarize, we find two channels for scalar instabilities. The first manifests as the presence of black hole hair, 
while the second (when the SLQL dimension becomes complex) corresponds to the bulk charged scalar condensing 
near the horizon of a black hole via Bose-Einstein condensation. 

For a spinor, a black hole hair gives rise to an isolated boundary Fermi surface, while when the SLQL dimension 
becomes complex, one finds a bulk Fermi surface in the near horizon region and the system settles into a Fermi liquid 
at very low energies. 

The differences between the scalar and spinor do not depend on any specific details and only on their statistics. 

F. SLQL as a universal intermediate-energy phase 

As discussed after (2.15), the charged black hole (2.14) provides a "universal" geometric description for many 
different theories at a chemical potential. Thus the SLQL phase appears as a universal IR fixed point among a large 
class of field theories at finite density, independent of their microscopic details. However, from equation (2.28), the 
extremal charged black hole has a finite zero-temperature entropy density. This implies that the ground states are 
highly degenerate. In a system without supersymmetry (as in this case 25 ) or other apparent symmetries to protect 
such degeneracies, we expect this nonzero entropy density has to with the large N limit we are working with. In 



Independent of the supersymmetry of the UV conformal field theory, the presence of the chemical potential is expected to break it. 
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other words, this nonzero zero-temperature entropy density likely reflects the existence of a large number of closely 
spaced states which are separated from the genuine ground state by spacings which go to zero in the N — > oo limit. 
One expects that the system should pick a unique ground state, which may or may not be visible in the large N 
limit. Given that there are many nearly degenerate low energy states, the precise physical nature of the ground state 
should be sensitive to the specific dynamics of the individual system. In a system with supersymmetry, then ground 
state degeneracy is indeed possible, and the above argument does not apply. There are many such examples in string 
theory. 

There is also another independent argument [92] based on density of states, which appears to indicate that any 
sort of scaling symmetry in the time direction only cannot persist to arbitrarily low energies. As we discuss in detail 
in Appendix E, however, there are various subtleties with the argument, which makes it hard to draw a general 
conclusion from it. 

Depending on the spectrum of charged or neutral matter fields (their charges and masses etc.), it is known that an 
extremal charged black hole suffers from various bosonic and fermionic instabilities [45, 88, 93-96], which will lead 
to a different lower energy state. Some of the instabilities have already been mentioned earlier and we will discuss 
them in more details in next two sections. The Einstein-Maxwell system (2.13) can also be generalized to an Einstein- 
Maxwell-dilaton system in which the charged black hole is no longer a solution. Instead, one finds a solution with 
zero entropy density and Lifshitz-type scaling in the interior (at T = 0) [97-100]. In all of these situations, there exist 
parameter ranges in which a charged black hole provides a good description for some region of the bulk geometry. 

It then appears that among a large class of systems with different microscopies and different lower energy states, 
the charged black hole appears as an intermediate energy state. Thus the SLQL phase may be considered a universal 
intermediate energy phase [44] which connects microscopic interactions with macroscopic, low energy physics, as 
indicated in Fig. 9. 26 It is characterized by an SX(2, R) symmetry including scaling in the time direction, but a finite 
spatial correlation length, as well as a nonzero entropy density. 

We should emphasize that for a given holographic system there may not always exist an energy range over the 
SLQL manifests itself as an intermediate state. SLQL behavior is manifest when there exists a hierarchy between 
the chemical potential \i and the energy scale at which a more stable lower energy phase emerges. On the gravity 
side there is then an intermediate region of the bulk spacetime which resembles that of AdS2 x (or its finite 

temperature generalization). In a situation where such a hierarchy does not exist, as we will discuss later, the SLQL 
can nevertheless provides a useful description for understanding the emergence of the lower energy state. 



FIG. 9. The semi-local quantum liquid phase is a useful description at intermediate scales; many different UV theories (i.e. 
those with gravity duals described by a Maxwell-Einstein sector) flow to it, and at low energies it settles into one of many 
different ground states, e.g. a Fermi liquid, a superconducting state, or an antiferromagnetic state, as will be discussed in 
subsequent sections. 

As emphasized by Anderson some time ago [102], the existence of a "universal intermediate phase" appears to be a 
generic phenomenon in nature; the familiar examples include liquid phases of ordinary matter, through which materials 
settle into different crystal structures at low temperatures. It is thus tempting to speculate that the appearance of a 





AFM 



SC 



The idea that an extremal black hole should be interpreted as an intermediate-energy state has been expressed by many people including 
the references [40, 42, 77, 78, 101]. 



29 



SLQL phase in holographic systems at a finite chemical potential may not be specific to these systems and may have 
wider applications. We will discuss more of these possible applications in the conclusion section. 

In the next two sections we will see examples of the evolution of the SLQL to various lower energy phases including 
Fermi liquids, superconductors, and AFM-like states. We find that the emergence of these lower energy phases can 
be characterized as a consequence of bound states formation in the SLQL. Conversely, SLQL can be considered as a 
universal deconfined (fractionalized) phase of these lower energy phases. 



III. HOLOGRAPHIC (NON-)FERMI LIQUIDS 

In this section we consider more explicitly the case when the holographic system has a Fermi surface and study the 
properties of small excitation around the Fermi surface. We will see that with a choice of parameters, the IR fixed 
point SLQL discussed in Sec. II C gives rise to non- Fermi liquid behavior which is the same as that of the strange 
metal phase of cuprate superconductors discussed in Sec. I A. However, when N is finite, we find that for most of the 
parameter range, the SLQL does not persist to arbitrarily small energies. Instead, at sufficiently small energies, the 
SLQL instead goes over to a Fermi liquid of heavy fermions through formation of fermionic bound states, where the 
scale at which this happens is exponentially small in N. 



A. Setup 

Consider a bulk Dirac spinor field ip(z,t,x) dual to a boundary theory spinor operator (D(t,x). For definiteness we 
will now specify to d = 3. We will assume a standard Dirac action in the bulk 



(Px^f^ i(t/jT M V M ^ - ra^), (3.1) 

where the covariant derivative T>m contains couplings to both the background gauge field and the spin connection of 
the charged black hole geometry (2.14), 

V M = d M + \u) ahM Y ah - iqA M . (3.2) 

Here q is the charge of O under the U{1) current J M . The asymptotic structure of solutions to (3.1) is different from 
the scalar case discussed above, essentially because it is a first order system with more components. In particular, 
the number of components of the boundary operator O is half of that of the bulk spinor [83], i.e. while ijj has four 
components, O has only two components, as expected for a three-dimensional field theory. The scaling dimension A 
of O is related to the bulk mass m of ip by 

A=^+mR. (3.3) 

We refer readers to [38, 40, 83] for detailed analysis of the Dirac equation and only mention the main results here. 
The boundary theory Green's function can be diagonalized to take the form 

G R (u;,k)=( Gl{uJ ' k) ° V (3.4) 
V ' \ G 2 (uJ,k) J V ' 

where the momentum k can be taken to be along one of the spatial directions and both G a =i y 2 are the ratio of 
appropriately defined boundary theory expansion coefficients A ai B ai as in (2.11) 

G a = ^. (3.5) 

The master formula (2.56) applies to both of the eigenvalues Gi 5 2, which will have different values for the prefactor in 
the AdS2 Green's function and the UV expansion coefficients a±,b±. A useful relation is Gi(cj, k) = G^cj, — fc); this 
means that with no loss of generality we can simply restrict to G2; when we refer to Gr in the rest of this section 
this is what we mean. For convenience we reproduce here (2.56), 



M^*0 + b_(uj,k)g k (uj)ii- 2 ^ 

a+(u), k) + a_(u;, k)Gk(u)li* 



Gnio;, k) = u T,i. rZ k T[ - 2 ,» . (3-6) 
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where Gk(w) is given by (2.52)-(2.53) and all of the UV coefficients a±,b± in (3.6) have analytic expansion in uj\ 

a+(u;,fc) = a { °\k) + a { +\k)u + a™ (k)u 2 + . . . . (3.7) 

When mR G [0, |) there is also an alternative quantization for the bulk spinor -0 which leads to the conformal 
dimension of the boundary operator O to be 

A = ^ — mR . (3.8) 

Below we will treat both quantizations together by allowing mR to take a negative value in (3.3). 

For simplicity below we will also set ji* to unity to avoid cluttering formulas. It can be easily restored on dimensional 
ground. 



B. Existence of Fermi surfaces 



As we discussed in Sec. HE 1, a boundary system has Fermi surfaces for a spinor operator 0, if there exists some 
fcp> ^\kp) = 0, i.e. the charged black hole admits an inhomogeneous fermionic hair. To see whether such a kp 
indeed exists one needs to solve the static Dirac equation (i.e. at uj = 0) for bulk spinor field ip in the full charged 
black hole geometry, which can only be done numerically. Squaring the static Dirac equation and rewriting it as a 
Schrodinger equation one can then reduce the problem of finding a Fermi surface to finding a bound state in the 
resulting potential as a function of k. It can then be checked explicitly that such a kp does exist for certain range 
of parameters (g, m) and when q is sufficiently large there can also be multiple bound states which corresponds to 
the system having multiple Fermi surfaces. In Fig. 10 we show the values of kp found numerically for a few choices 
of mR, taken from [40]. In the figure, the oscillatory region k < k Q (see (2.81)), where becomes pure imaginary 
and Gr becomes oscillatory in \oguj (see (2.82)), is shaded. As discussed in Sec. HE 3, the oscillatory region can be 
attributed to a bulk Fermi surface in the near horizon AdS2 region, which gives rise to a "Fermi ball" in the boundary. 
In contrast, the Fermi surface at kp from a fermionic hair is isolated. 

For standard quantization, i.e. m > 0, the parameter region in which a kp exists always lie inside the parameter 
region (2.81) for which the oscillatory region exists. But note that kp is always greater than fc , i.e. lying outside the 
oscillatory region in momentum space. For alternative quantization, there exists a parameter range in which we have 
a single isolated kp even in the absence of an oscillatory region. 




FIG. 10. The values of kF as a function of q for the Green function G2 are shown by solid lines for mR = 0, 0.4. In this plot 
and ones below we use units where R= 1, /i* = 1, #f = 1 and d — 3. The oscillatory region, where Vk — ^ \jk 2 + m? — ^ is 
imaginary, is shaded. Given G\(k) = G2(—k), so kF for G\ can be read from these plots by reflection through the vertical k — 
axis. For convenience we have included in each plot the values of kF for the alternative quantization using the dotted lines. 
Thus the dotted lines for the second plot (mR = 0.4) correspond to mR = —0.4. Also for m = the alternative quantization is 
equivalent to the original one. This is reflected in the first plot in the fact that the dotted lines and solid lines are completely 
symmetric. Both plots are symmetric with respect to q, k — > —q, —k. 
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FIG. 11. Distribution of v^ F in the m — q plane. When there are multiple Fermi surfaces, v^ F is calculated for the one with 
the largest k F - The two plots correspond to G2 and G\ components. In the white region, there is no Fermi surface. There 
is no long-lived quasiparticle in the orange region, Vk F < |. In the remainder of the parameter space, there are log-lived 
quasiparticles. Above the dashed lines there is an oscillatory region. 

C. Fermi surfaces with or without long-lived quasiparticles 

Near a k F and at low frequencies, the retarded function Gr is given by (2.66)-(2.68), which we also reproduce here 
for convenience 



G R (uj,k) 



with 



v F 



k — k F 



d k a+\k F ) 
J 1 ) 



-u ■ 



((M 



hi 



E(w) = hg kF (oj) = hc(v kF )u? l ' k i 



hf(k F ) 



h 



(k F ) 



d k a^\k F ) 



(3.9) 



(3.10) 



d k a^\k F y 

The self-energy E depends on frequency only and is characterized by a nontrivial scaling exponent v kp , which can be 
computed from (2.51) once k F determined. The result is presented in fig. 11. 
Poles of (3.9) in the complex cj-plane, which are determined by 

k-k F - —u c - E(cj c ) = uo c {k) = uojk) - iT(k) (3.11) 

v F 

give the dispersion of small excitations near the Fermi surface. The real part cj*(fc) of u c (k) gives the dispersion 
relation and the imaginary part T gives the width of an excitation. Given that in (3.11) only E is complex, the width 
T is controlled by the imaginary part of E which is in turn is proportional to ImQ kF (oo). The physics turns out to 
be rather different depending on the value of v kp , and we now study the different values in turn. For simplicity of 
notation in the rest of this section we will drop the subscript k F : v = v kp in this section only. 



1. v>\ 



If v > ^, then the analytic linear-in-cj term in (3.9) dominates over the self-energy, which goes like uj 2u . Thus the 
dispersion is essentially controlled by the balancing between the first two terms in (3.11); note however that the sole 
contribution to the imaginary part of the pole comes from the complex self-energy. We find the following expression 
for the motion of the pole as k is varied near k F : 

2v 



u c (k) = v F (k — k F ) — v F hc(v) 



v F (k — k F ) 



(3.12) 



Thus we see that this has a linear dispersion relation uo* = v F {k — k F ), with width T ~ uj 2v that is always much 
smaller than the energy uj*: as we approach the pole ~ (k — /c^) 2 ^ -1 as - k F —^0. Furthermore the 

residue at the pole is always finite and given by 

Z = -h t v F . (3.13) 

This is then a long-lived quasiparticle. It is not qualitatively different from the excitations of a Fermi liquid, except 
that its width satisfies T ~ uj 2v rather than ~ uj 2 . 
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2. v < \ 

If v < ^, the situation is different. Now in (3.9) the non-analytic contribution from E(u;) always dominates over 
the analytic linear-in-cj correction, which we may completely neglect. The solution to (3.11) now takes the form 

fk-k F \^ , AS 

u>c = t**[ -TTT- ( 3 - 14 ) 



hc(v) 

As c(y) is complex this pole has both a real and imaginary part; however now they are both coming from the same 
place, i.e. the AdS2 Green's function. As a result the dispersion is now nonlinear, uu* ~ (k — kp)^ . Importantly, the 
width T is now always comparable to the frequency cj*, e.g. for k — kp > we have 

r(fc) /arg(c(i/))\ /0 „ r , 

- tan 6V v y; = const (3.15) 



tJ*(fc) V 2l/ 

A similar expression exists with a different angle for k — kp < 0. Thus as we vary k — kp through 0, the pole moves 
along a straight line towards the origin, eventually bouncing off and retreating at another angle (but always in the 
lower-half plane). The residue at the pole is 

UcHl (k-k F )^~\ (3.16) 



2v(k - k F ) 

and so actually vanishes as we approach the Fermi surface k — kp = 0. The fact that the width of the excitation is 
always the same order as its energy and that the residue vanishes at the pole indicate that this can not considered as a 
quasiparticle; thus we see that we have an example of a system with a sharp Fermi surface, but with no quasiparticles! 

3. v — - : Marginal Fermi Liquid 

We turn now to the threshold case between the two studied above. Here we see that the two terms (analytic and 
AdS2) are both naively proportional to uo and are thus becoming degenerate. As is often the case when there is such 
a degeneracy, a logarithmic term in uj is generated. Basically in (3.9), v F has a simple zero and c(v) in the self-energy 
E (see (2.52)) has a simple pole at v = \\ the two divergences cancel and leave behind a finite uloguj term with a 
real coefficient, 

G R (w, k) « - h \ (3.17) 

k — hp + ciu log uo + Ciuj 

Here c\ is real and c\ is complex; for more details on the derivation of this formula see Section VI A of [40] . Remarkably, 
this is of precisely the form postulated in [12] for the response function of the "Marginal Fermi Liquid" discussed earlier 
around (1.13). Here this response function emerges from a first-principles holographic calculation (albeit only at the 
specially tuned value v = |). Note that here the width is suppressed compared to its energy, but only logarithmically 
r ~ Similarly, the residue of the pole scales like Z ~ 1/ log(fc — kp) and so also vanishes logarithmically. The 

potential physical applications of this v = | point make it very interesting, although from our gravitational treatment 
so far we do not really have any way to single it out. 

4- Perspective from low energy effective theory 

The qualitative behavior discussed above, in particular, the presence (absence) of long-lived quasiparticles for v > \ 
(y < |) can also be seen from the low energy effective theory (2.75) based on simple dimensional analysis [43]. Since 
in the SLQL, k does not scale (i.e. has dimension 0), from (2.76) \£ should have scaling dimension in the IR. Thus 
the hybrization term in (2.75) has dimension | + v (which is the dimension of <£). As a result the hybridization 
coupling is relevant for v < | and irrelevant for v > |. When the hybridization coupling is relevant, free fermion ^ is 
strongly mixed with <I>, resulting a large decay width and a breakdown of the quasiparticle picture. When v — \, the 
coupling is marginal. Note that while the logarithmic suppression in the decay width was the original motivation for 
the word "Marginal" in the "Marginal Fermi Liquid" given in [12]; here we find it is due to a marginal hybridization 
coupling and the word "Marginal" is entirely appropriate. 
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D. Transport of holographic Fermi surfaces 

Having understood the single-particle response functions of our system, we now turn to collective properties, such 
as transport. A very natural observable to compute is the conductivity; in particular, for Landau Fermi liquids this 
is essentially always controlled by the single-particle lifetime, and we find for the DC conductivity <jfl ~ T 2 . What 
is <t(T) in our model? 




O(N ) 



+ + + + + + + + 



FIG. 12. Conductivity from gravity. The current-current correlator in (3.18) can be obtained from the propagator of the gauge 
field A x with endpoints on the boundary. Wavy lines correspond to gauge field propagators and the dark line denotes the 
bulk propagator for the tp field. The left-diagram is the tree level 0(N 2 ) contribution while the right diagram includes the 
contribution from a loop of ip quanta. The contribution from the Fermi surface associated with boundary fermionic operator 
O can be extracted from the diagram on the right. 

To understand the calculation to follow, it is necessary to understand the answer to the following question: how 
much of the charge of our system is carried by the Fermi surface that we have just identified? Note from (2.20) 
that the background charge density of our system scales like p ~ ^2, where ^ is the inverse Newton's constant 
that multiplies the whole bulk action. In known examples of the duality, this is typically proportional to a positive 
power of N, i.e. N 2 or N?; for simplicity let us just assume the former and say that p ~ 0(N 2 ). But k F identified 
above has no factors of n 2 in its definition; thus we conclude that kp ~ 0(N°). Since this is the Fermi surface for 
a singlet operator, it means that the charge density that is associated with our Fermi surface is also 0(N°). Thus 
we conclude that our Fermi surface makes up a tiny (i.e. large-N -suppressed) portion of the full charge density. In 
our gravitational description, the rest of the charge is stored behind the black hole horizon; in field theory terms, this 
presumably means that it is stored in some other degrees of freedom about which we know very little except that in 
some sense there are "TV 2 of them" . 

This means that none of the leading thermodynamic or transport properties will be sensitive to the Fermi surface. 
We can directly see this from a bulk calculation. The conductivity can be found from the Kubo formula, 

a(u) = ^{J x {bj)J y {-Uj)) retarded- (3.18) 

We can find this retarded response function from gravity using the procedures described above; however if we do that 
it is clear that the leading answer will be of 0(N 2 ) and will depend only on classical fluctuations of the gauge field, 
which know only about the black hole and appear to know nothing about the Fermi surface. It is clear that to see 
the effects of the Fermi surfaces identified above, we have to go to higher order in 1/iV 2 ; this maps to quantum effects 
in the bulk. Thus we must perform a one-loop calculation on the gravity side, including a fermion loop in the gauge 
boson propagator as shown in Figure 12. 



1. DC conductivity 

In many ways this calculation parallels that of Fermi liquid theory, except for the extra bulk dimension and com- 
plications from the curved geometry. A main difference comes from the fact that we are performing the computation 
in a black hole spacetime, meaning that processes such as Figure 13 where the fermion goes through the horizon can 
contribute. We do not provide any details of the computation here, only pointing out that the final result can be 
written in terms of boundary theory quantities as 

ct(uj) = - fdk [ ^1^2 /K)-/M A(u)i,k)A(cji,uj2,v, k)A(oj 2 ,ui,w, k)A(uj 2 , k) (3.19) 

ZUJ J J Z7T Z7T OJi — UJ — Co>2 — ie 
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FIG. 13. Possible contribution to conductivity in which fermion loop passes through the horizon. This represents the decay of 
the current due to current-carrying particles decaying by falling into the black hole, and is the dominant contribution to the 
DC conductivity. 

Here A = ^IirGr is the single-particle spectral function described in the last section, f(u) = is the Fermi 

distribution function. All radial integrals have been absorbed into an effective vertex A, which takes the schematic 
form A ~ J dr t/j(r)KA(r)^(r)^ where ip(r) is a fermion wavefunction and Ka{?) is a bulk-boundary propagator for 
the gauge field (which also includes a contribution from the graviton, which mixes with it on this background). This 
expression can be interpreted in boundary theory terms, as shown in Figure 14; the effect of the radial direction is to 
provide an exact expression for the renormalized vertex A. 



A((i),k) 




A(-(»,k) 

FIG. 14. Interpretation of (3.19) in field theory terms; the conductivity can be written as an expression similar to that in Fermi 
liquid theory, but with exact propagators and renormalized vertexes that can be found from gravity. 

In general A is a complicated function of all its arguments, but in the low temperature limit and near the Fermi 
surface it becomes a smooth real function of |fc|, independent of uj and T. In this case (3.19) shows that the conductivity 
is controlled by the single-particle spectral function near the Fermi surface, and we find the expression 

a DC = a(q,m)T- 2vk F. (3.20) 

Here a(q,m) is a non-universal number that depends on a radial integral over the bulk fermion wavefunctions (and 
can be found numerically) but the exponent in T is determined by the single-particle decay rate and so by the SLQL. 
In particular, note that for the Marginal Fermi Liquid case v\ ZF = ^, we find that the contribution to the resistivity 
from the Fermi surface is linear, p ~ T, as observed for the cuprates in the strange metal phase and many other 
non- Fermi liquid materials [10]. 

Note that in general the electrical transport scattering rate - which emphasizes large momentum transfer - can be 
different from the single-particle scattering rate. On the other hand in our system the exponent in (3.20) is determined 
by the single-particle scattering rate. On the gravity side this arises from the fact that the dissipative portion of the 
correlation function is determined from Figure 13 essentially by how fast the bulk fermions fall into the black hole, 
which determines the single-particle decay rate and knows nothing of momentum. 27 On the field theory side this is 
due to the fact that the SLQL sees only the time direction and depends very weakly on momentum. 

2. Optical conductivity 

In the computation above the formula (3.20) for the DC conductivity is rather insensitive to whether or not the 
Fermi surface has long-lived quasiparticles or not , i.e. the value of v^ F . The optical conductivity <j{oo) can also be 



Note that higher loop diagrams which could potentially spoil this equivalence are suppressed by higher powers of 1/N 2 . 
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computed along similar lines; this turns out to depend critically on the value of v^ F . Again, we present here only 
results. 

For Vk F > \ there are two regimes. In the first we hold ujT~ 2uk F fixed in the low temperature limit and find then 
that we can approximate <j(uj) by a Drude form, 



uj 2 



a(uj) = - p , r- 1 = 2v F ImS(cj = 0) - T 2vk ? (3.21) 



• iuj 



with uj 2 an overall constant. This is in fact the standard answer for a weakly coupled theory with well-defined 
quasiparticles, whose lifetime is given by r. This should not be surprising; in some sense the effect of holography (in 
this regime) is to provide the nontrivial lifetime r by coupling these quasiparticles to a quantum critical bath. In the 
other regime uj ^> T (but still uj <C /i) we find 



uj 2 



(j(uj)=i^+b{iuj) 2vk F- 2 , (3.22) 



UJ 



with b a real constant. The 1/uj term gives us a term proportional to 5(uj) with weight consistent with (3.21). 
On the other hand, for v^ F < |, when we have no long-lived quasiparticles, we find a scale-invariant form, 

aH = T" 2 ^F 1 (^), (3.23) 

where F\(x) is a scaling function that is constant as x — >> 0. For large x it falls off as x~ 2vk F , giving us an optical 
conductivity a(uo) ~ a(iuj)~ 2vk F with a a real constant. This behavior is consistent with a system with no scales and 
no quasiparticles. 

Finally, for the marginal case u^ F = |, we find behavior similar to (3.23) but with logarithmic violations, i.e. 

a{uj) = T- 1 F 2 (jpl°g^) (3-24) 
where again F2 goes to a constant at small uj. For uj ^> T in this case we find 



E. SLQL as a fractionalized phase of Fermi liquids 

In this subsection we further explore the physics of the oscillatory region k < k Q discussed in Sec. HE 3. Such a 
region exists when 

^f- > m 2 R 2 (3.26) 

and it is shown in the shaded regions in Fig. 10. In the oscillatory region, the IR scaling dimension of a fermionic 
operator is complex and the spectral function is nonzero at uj = 0. At a nonzero uj, the correlation function (2.82) is 
periodic in \oguu and has an infinite number of poles in the lower half cj-plane. It can also be readily checked that all 
the poles have comparable real and imaginary parts and thus cannot be considered as corresponding to quasiparticles. 
From the bulk point of view, this region can be interpreted as corresponding to the projection of a bulk Fermi surface 
(which is one dimension higher) to the boundary. Small excitations at the bulk Fermi surface have a large decay rate 
since they can fall into the black hole without experiencing any potential barrier. 

We will now show that in the parameter range of (3.26), at any finite N the system in fact goes over to a Fermi 
liquid of "heavy" fermions at low energies. We will keep the boundary spacetime dimension d general. 

From (2.89), the local charge density of fermions is given by 

PL=D d {k ( P)\ D d = qs^y d (3.27) 

where is the volume of a unit sphere in d-dimension and s is the number of possible spin degeneracies. pL in (3.27) is 
finite and scales as 0(N°) in the large N limit 28 which naively implies that the fermionic backreaction can be ignored 



In contrast the charge density carried by the black hole, given in (2.28) scales as 0(N 2 ). 
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at leading order. However, the contribution to the boundary charge density from pL is obtained from integrating 
it over the radial direction of the black hole, and since in the near horizon region (2.30) the proper radial distance 
is infinite and the local proper volume of the transverse R d_1 is constant, the total boundary density coming from 
integrating (3.27) over the radial direction is in fact infinite [95]. 



AdS 2 x R 2 AdS 4 




Lif z AdS 4 



FIG. 15. Two different geometries (with d = 3): on the top, the AdS2 x R 2 describing the SLQL phase; on the bottom, its 
resolution into a Lifshitz geometry with a finite z given by (3.29). The horizon direction represents the y direction, while the 
vertical direction represents the transverse R 2 . In the plot for the Lifshitz geometry it should be understood that the tip lies 
at an infinite proper distance away. When z is large as in (3.29), there is a large range of y for which the Lifshitz geometry 
resembles that of AdS2 x R 2 . Also note that e~ y translates into the boundary theory energy scale. In the Lifshitz geometry, 
the black hole has disappeared; it is replaced by a naked singularity with zero horizon size and zero charge. All the charge 
density is now carried by a bulk fermionic gas and the entropy density vanishes at zero temperature. 

The backreaction of the fermionic gas is thus important no matter how large N is. The backreacted geometry can 
be found by using (3.27) and its corresponding stress tensor as the source for bulk Maxwell and Einstein equations. 
We refer readers to [95] for details of the calculation and only state the results here. One finds that after taking into 
account of the backreaction the near-horizon geometry is modified from (2.87) to a Lifshitz geometry 

ds 2 = -e-^dt 2 + e~^^R 2 dx 2 + R 2 dy 2 , A t = ^e"* (3.28) 

with 

where we have only kept leading order terms in 1/z expansion. It is important to emphasize that even in the large 
N limit, the exponential prefector before the second term of (3.28) cannot be set to unity as it becomes important 
when y is sufficiently large (i.e. y ~ O(z)). Also note that as (2.87), (3.28) should be used in the region y > 0; the 
geometry matches to the rest of the rest of black hole (UV region) near y ~ 0. 

In the backreacted geometry (3.28), shown schematically in Figure 15, the local proper volume of the transverse 
R d_1 goes to zero as y — >• 00, resulting in a finite boundary fermionic density of order 0(7V 2 ), given by 

Pf = R2(P*R) PL J dye * y = . (3.30) 



Plugging the explicit value (3.29) of z into (3.30) we find that (3.30) becomes identical to the total charge density (2.28), 
i.e. all the charge density of the system is now carried by the fermioinic gas. The black hole has disappeared! It is 
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now replaced by a naked singularity 29 of zero horizon size and zero charge. The system now also has zero entropy 
density. Note that such a singularity appears to be of the "good" type allowed by string theory [103, 104]. As we will 
see below, things could still fall into a singularity, which means in the boundary it should be interpreted as some of 
"bath" which can dissipate things away. While it has zero entropy and charge at order 0(N 2 ), likely such a "bath" 
still carries some charge and entropy at subleading order. 

Recall that in obtaining (2.89) and (3.27) we used the approximation of a local chemical potential (i.e. Thomas- 
Fermi approximation), which is valid when qgF-, mR, R are taken to be parametrically large. Thus (3.29) and (3.30) 
are also only valid in this regime. 

The Lifshitz geometry (3.28) is well approximated by (2.87) until y ~ z. From the boundary theory perspective, 
this implies that at an energy scale of order e~ z ~ e~°( N \ physics will start deviating from that described by the 
SLQL. In other words, at sufficiently low energies, the SLQL will order into a new phase which is described by the 
Lifshitz geometry. What is this new phase? It turns out to be a Fermi liquid of "heavy" fermions [44]! Below we 
outline the basic arguments, following the treatment of [44], in which more details can be found. Similar calculations 
appear in [107, 108]. 
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FIG. 16. Plots of the WKB potential (3.31) for various values of uj in the Lifshitz region. For uj = 0, the potential is bounded 
from below and the Bohr-Sommerfeld quantization gives normalizable modes in the bulk which correspond to different Fermi 
surfaces in the boundary theory. For an exponentially small uj, although the potential is unbounded from below in the large y 
region, the excitations have a small imaginary part due to the potential barrier. For small uj but not exponentially small, the 
potential barrier disappears and the potential becomes the same as that for the AdS2 x R 2 . 

The basic idea is again to examine the retarded Green function for the fermionic operator in which we now expect 
to find new physics at sufficiently low energies. For this purpose one needs to solve the Dirac equation for the bulk 
fermions in the backreacted geometry (3.28), which for q and mR parametrically large can be reduced to a Schrodinger 
problem describing a particle of zero energy moving in a potential 

V(y) = -L (k 2 e^ -k 2 Q - uj 2 e 2 y - y/l^que*) . (3.31) 

For the AdS2 x geometry, the first term is simply k 2 . 

The final results can be readily understood from the qualitative features of the WKB potential (3.31), as shown 
in Fig. 16 for k 2 < k 2 . First consider uj = 0, the potential starts being an almost negative constant for y <C z, but 
eventually rises above zero and approaches +oo as y — > +oo. The Bohr-Sommerfeld quantization condition then leads 
a series of normalizable bound states at discrete values of k. Each such normalizable mode then gives rise a boundary 
Fermi surface. Alternatively, from the bulk point of view, we now have a bounded potential governing movement 
in the radial direction. This quantizes the radial modes, and the continuous bulk Fermi surface indicated in Fig. 8 
becomes discretized. When projected to the boundary this now become a series of Fermi surfaces as indicated in 
Fig. 17. Now notice that in (3.31), for any nonzero uj no matter how small, the a;-dependent terms will eventually 



Note that at this "singularity" all scalar curvature invariants are actually finite, but nevertheless an observer freely falling into the 
singularity will feel divergent tidal forces [105]. 
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dominate for sufficiently large y and drive the potential to — oo as y —> +00. This implies bulk particles will eventually 
fall into the singularity of the Lifshitz geometry and the corresponding boundary fermionic excitations should develop 
a nonzero width. For exponentially small a;, there is still a potential barrier and the tunneling rate is exponentially 
small. In this range there are long-lived quasiparticles. More explicitly, one finds that near each kp the retarded 
function for uj in the range 

k ( k \ 

uj <C oj c (k) = — exp I — z log — - I (3.32) 



z \ kp 

can be written as 



with 



G R (u>, k) = (3.33) 

uj — v F (k — k F ) + 2J 

Z, vp ~ exp ( — z log , £ ~ i exp | — ^ Z ^ F (hi ^\ — z \og^-\ . (3.34) 

V k F J \ M \u J k F J 

Note that since z ~ 0(7V 2 ), these are very heavy fermions 30 and the decay rate is exponentially small in uj. 

For small but not exponentially small a;, the first term in (3.31) is not important for any value of y and the potential 
barrier disappears. One thus recovers previous results of SLQL. In particular, the decay rate becomes O(l) and there 
are no long-lived quasiparticles. 

This procedure of solving the Dirac equation in the backreacted geometry (3.28) is self-consistent and extracts 
the leading non-perturbative behavior in 1/N 2 . There are also perturbative loop corrections in the bulk, which give 
perturbative corrections in 1/N 2 to the self-energy. We expect the qualitative features of the above discussion (e.g the 
family of densely spaced Fermi surfaces etc) to be robust against these corrections, as they have to do with the global 
structure of the backreacted geometry (3.28). Near the Fermi surface perturbative corrections to the self-energy will 
give rise to a term jf^uj 2 with c some complex 0(1) coefficient, which will dominate over the imaginary part of E 
in (3.33) for uj in the range (3.32). Thus we expect that the quasiparticle decay rate should be proportional to uj 2 as 
in a Landau Fermi liquid. 31 

We thus see that at low energies the oscillatory region k < k Q splits into a large number of closely spaced Fermi 
surfaces whose excitations are Landau quasiparticles. In this regime the Luttinger theorem should hold, which can 
be checked explicitly (see [44] for details). 

Different Fermi surfaces in the family of densely spaced Fermi surfaces can be interpreted as corresponding to 
different bound states generated by the fermionic operator 0, as each of them corresponds to a different radial mode 
in the bulk. Recall that in our set-up, O is a composite operator of fundamental fields. The number of degrees of 
freedom for the fundamental fields is 0(N 2 ). Thus the Fermi liquid state can be considered a "confined" state, in 
which the low energy degrees of freedom are Fermi surfaces from a discrete set of composite fermionic bound states. In 
contrast, the SLQL is a "deconfined" state in which the composite bound states deconfine and fractionalize into more 
fundamental degrees of freedom. Presumably these more fundamental degrees of freedom account for the entropy 
density of the SLQL. In the bulk, the emergence of the fractionalized SLQL phase is reflected in the emergence of a 
charged black hole description. In the Fermi liquid state, the system is characterized by "heavy" fermions. But such 
coherent quasiparticles disappea in the SLQL. Instead one finds some kind of quantum soup which is characterized 
by scaling behavior in the time direction for any bosonic or fermionic operators. In the current context the presence 
of a large number of Fermi surfaces has to do with a spectrum of densely spaced bound states. The fractionalized 
picture should be independent of this feature. We note also that this interpretation of the SLQL as a fractionalized 
phase resonates with the lattice model construction of [77], although the details are different. 

Note that the WKB potential (3.31) only includes the near horizon region of the backreacted charged black hole 
geometry, and does not include the Fermi surfaces discussed in Sec. Ill B. In Fig. 18 we show a cartoon of the WKB 
potential for uj = for the full spacetime, i.e. including the asymptotic AdS^+i region. The isolated Fermi surfaces 
discussed earlier in Sec. IIIB and IIIC appear as bound states in a potential well outside the near-horizon region. 
Such isolated Fermi surfaces exist both in the Fermi liquid phase and in the SLQL, i.e. the corresponding fermionic 
excitations remain confined in the deconfined SLQL phase. At exponentially small energies, the excitations again 
become Landau quasiparticles, but with an O(l) effective mass. 



The exponentially small Fermi velocity vp implies that the effective mass ra* is exponentially large, i.e. ra* ~ exp y-\-z\og 
Nevertheless, one should keep in mind that the much smaller non-perturbative correction does signal some nontrivial underlying physics 
beyond that of a Landau Fermi liquid. 
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FIG. 17. In the Lifshitz geometry, we find a set of discrete states in the bulk in which radial motion is quantized. This results 
in a family of concentric Fermi surfaces in the boundary theory, which resolves the Fermi disk of Fig. 8. 
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FIG. 18. A cartoon of the bound-state structure for the full spacetime, i.e. including the asymptotic AdS4 region. The Fermi 
surfaces in Sec. Ill B and III C appear as bound states in a potential well in the UV region. Note that the plot only intends to 
indicate the locations of bound states in the bulk spacetime and does not reflect the genuine WKB potential. 

To summarize, imagine a system with only a single fermionic operator O satisfying the condition (3.26), and no 
other instabilities. Then at very low energies, the system is described by a Fermi liquid with 0(N 2 ) densely spaced 
Fermi surfaces, each of size 0(N°). The quasiparticle excitations have a very large effective mass (proportional to 
e N ). There could also be some additional isolated Fermi surfaces with an O(l) effective mass. At small but not 
exponentially small frequencies, there is a wide energy range over which the system is controlled by the SLQL, with 
scaling behavior in the time direction and various non-Fermi liquid behavior. In the SLQL phase the clue that the 
system will eventually settles into a Fermi liquid state is the existence of a region k < fc D , where the scaling dimensions 
for fermions become complex. 

F. Generalizations 

Before concluding this section, we mention various generalizations which we do not have space to cover. 

In the presence of a magnetic field, one can show that fermionic excitations fall into Landau levels [106, 109, 110] 
and that the free energy displays de-Haas van Alphen oscillations in the inverse magnetic field [111, 112]. The latter 
calculation also requires a one- loop calculation in the bulk. In [112] it was shown that while the oscillations exist 
with the expected period, for v^ F < | the amplitude of the oscillations does not take the textbook Lifshitz-Kosevich 
form. The fermion response in various holographic superconducting phases has been studied in [113-119]. The 
superconducting condensate opens a gap in the fermion spectrum around the chemical potential if a certain bulk 
coupling between the spinor and scalar is included [114]. In particular, in the condensed phase, one finds stable 
quasiparticles, even when v^ F < \. The properties of Fermi surfaces and associated excitations with more general 
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fermionic couplings or boundary conditions have also been studied in [120-123] and in more general geometries in 
[124, 125]. Fermi surfaces in 1 + 1-dimensions were studied in [126-128]. References [129, 130] considered the pairing 
instability of holographic non-Fermi liquids. 

By considering spinor fields whose mass and charge grows with N (i.e. m,q ~ N) and a gauge coupling of order 1/iV, 
the Lifshitz geometry (3.28) has been generalized to have a z ~ 0(N°) in [131] (see also [132-134]). Various aspects 
of the boundary system dual to such an "electron star" were studied in [107, 108, 135-137]. In particular, it was 
found [107, 108] that the fermionic spectral function exhibits multiple Fermi surfaces with exponentially suppressed 
(in frequency) decay rates very similar to the behavior discussed earlier around (3.32)-(3.34) (see also [138] for a 
hard- wall model for constructing holographic Fermi liquid). 



IV. SCALAR INSTABILITIES AND QUANTUM PHASE TRANSITIONS 

The previous section discussed fermions in some detail, describing both non-Fermi liquid physics and the ultimate 
Fermi- liquid- like ground state of the system. In this section we turn to bosons. Recall from Section HE that while 
fermionic response functions can have gapless modes that describe propagating degrees of freedom (i.e. are in the lower 
half complex -uj plane), generically for bosons we instead find instabilities, i.e. poles in the upper half complex plane. 
This indicates that the scalar operator is trying to condense, and if the operator is charged under any symmetries 
then they will be broken in the condensed phase. 

In this section we examine these bosonic instabilities of the finite charge density system. In particular, when there 
is an instability we describe how to tune external parameters to make this instability vanish; if this can be done at 
zero temperature, then this defines a quantum critical point that separates two phases. We will study in detail the 
critical behavior near such a quantum critical point. We will mainly follow the discussion of [45-47] (related results 
have also been obtained in [139, 140, 142], see also [143]). For concreteness we will also specialize to d = 3. 

We first remind the reader of our setup; consider a CFT 3 with a scalar operator O. We discussed above in Section 
II D how to use holography to compute the retarded Green's function of O. Before moving on, we will first slightly 
generalize that discussion to also include double trace deformations in the dual CFT3 

SCFT 3 "> S C FT 3 + JO 2 . (4.1) 

It is straightforward to include the presence of such a deformation in our holographic computation. Details are 
provided in Appendix C; the main result can be found in (C6), from which we see that (2.56) becomes 



g r ^ %) = : + r r~ \T/\-2* k ( 4 - 2 ) 



, 2 ^ j+M + Mfc,^)g fc (^)/i- 2 ^ 

where 

a±(oj,k) = a±(uj,k) + K+b±(u,k) = (k) + 0(u) , a ( ±\k) = a { ± ] (k) + K+b ( ±\k) . (4.3) 

Note that in general the coefficient ft + of the double-trace operator provides a "knob" that can be turned to tune the 
system through a quantum phase transition. 

To follow the standard terminology for discussing phase transitions below we will start calling (4.2) the susceptibility 
and use the following notation: 

X = G r (oj = 0,fc = 0), X (k) = G R (oj = 0,k), X (^,k)=G R (uj,k) . (4.4) 

where we distinguish the three only by their arguments. From (4.2) we find the momentum-dependent and uniform 
static susceptibilities are given by 

and the full dynamical susceptibility takes the form 

X ( U , k) = & Mfc,") + Mfc,")foH^ m 

a + (k,uj) + a-{k,uS)Qk{w)lJi* k 

Below we first recall the instabilities discussed in Sec. HE for the system at a finite chemical potential and the 
corresponding quantum critical points. We will then discuss quantum critical behavior near these quantum critical 
points. 
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A. Finite density instabilities 



Recall that in HE we discussed two possible channels for scalar instabilities, which manifest themselves as presence 
of poles of (4.6) in the upper half complex-o;-plane: 

1. The SLQL scaling dimension \ + becomes complex for some /c, for which there are an infinite number of 
poles in the upper half u;-plane. Writing (2.38) as 




Vk = \ u + ^r, u = m 2 R\ + - - <£ (4.7) 



1 

4 

Vk becomes imaginary for sufficiently small k whenever u < 0. For a given m, this always occurs for a sufficiently 
large q. For a neutral operator q = 0, u can be negative for m 2 R 2 lying in the window 

--<m 2 R 2 <-l (4.8) 
4 2 v ) 

where the lower limit comes from the stability of vacuum theory (i.e. BF bound of AdS4) and the upper limit 
comes from the condition u < after using the relation (2.31). 

2. a+\k) can become zero for some special values of momentum kp. At k = kp it is clear from (4.2) that since 

ct+ = 0, x nas a singularity at uj = 0. Furthermore since at+ changes sign near k = kp, as discussed in 
Sec. HE 1, a pole moves from the upper half cj-plane (for k < hp) to the lower half cj-plane (for k > kp) through 
uj = 0. In Fig. 19 we show some examples of a neutral scalar field for which has a zero at some momentum. 




FIG. 19. Left: a ( +\k) and (k) plotted for different values of m 2 R 2 with q = 0; blue is m 2 = -1.4999, purple m 2 R 2 = -1, 

green is m 2 R 2 = —0.5, red is m 2 R 2 — and the black line is m 2 R 2 — 0.5. a+\k) is positive and monotonically increases with 

/c, while b+\k) monotonically decreases with k. Thus for ft+ < 0, a+\k) = a+\k) + K+b+\k) is a monotonically increasing 

function of k. Note that b+\k = 0) = for m 2 — 0. This implies that when m 2 > 0, a+\k) > for any ^+ < 0. Right: a+\k) 

for m 2 R 2 = -1.4999 and q = 0. (k) has a zero for some k when ft+ < k, c = —2.13. For > ft+ > k, c there is no instability 

(see ft+ = — 1 curve). When ^+ > 0, (k) can again develop a zero for some /cf (with the value of kF approaching infinity 
for ^+ — > 0+); this is a UV instability which is already present in the vacuum. 

In the parameter range (say for m, q, ft + ) where either (or both) instability appears, the system is unstable to the 
condensation of operator O (or in bulk language the condensation of </>). For a charged scalar the condensed phase is 
the well-studied holographic superconductor [93, 94]; for reviews see [51, 52]. As was first pointed out in [87], the first 
instability essentially drives the original construction of the holographic superconductor. Holographic superconductors 
due to the second type of instability were first discussed in [46] . 

For a neutral scalar, the first type of instability was first pointed out in [88], and as discussed in [45] the condensed 
phase can be used as a model for antiferromagnetism when the scalar operator is embedded as part of a triplet 
transforming under a global SU (2) symmetry corresponding to spin. For a single real scalar field with a Z2 symmetry, 
the condensed phase can be considered as a model for an Ising-nematic phase. Instabilities of the first type also 
appear in various other holographic constructions [139, 140]. The potential interaction between these two types of 
order is discussed in [141]. 
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Both types of instabilities can be cured by going to sufficiently high temperature; there exists a critical temperature 
T c , beyond which these instabilities no longer exist and at which the system undergoes a continuous superconducting 
(for a charged scalar) or antiferromagnetic (for a neutral scalar) phase transition. As have been discussed extensively 
in the literature [45, 89, 90] such finite temperature phase transitions are of mean field type due to that the boundary 
conditions at a finite-temperature black hole horizon is analytic. 

Alternatively one can continuously dial external parameters of the system at zero temperature to get rid of the 
instabilities. 32 The critical values of the parameters at which the instabilities disappear then correspond to quantum 
critical points (QCP) where quantum phase transitions into a superconducting or antiferromagnetic phase occur. We 
now quickly summarize the nature of each of these quantum critical points. 



B. Quantum critical points 

1. Bifurcating quantum critical point 

For the first type of instability a quantum critical point occurs when we violate the Breitenlohner-Freedman bound 
of the AdS 2 region at k = [45, 139], i.e. from (4.7), at 

u = u c = . (4.9) 

For example for a neutral scalar field (with q = 0) this happens at 

m 2 c R 2 = -\. (4.10) 

Note that while in AdS/CFT models the mass square m 2 for the vacuum theory is typically not an externally tunable 
parameter, the effective AdS 2 mass square can often be tuned. For example, if we place a holographic superconductor 
in a magnetic field at zero temperature the AdS 2 mass can be changed by tuning the magnetic field [45]. In the set-up 
of [139], a somewhat different order parameter is dual to a bulk field that has a mass that depends on the external 
magnetic field. See also [45] for a phenomenological model. Here we will not worry about the detailed mechanism to 
realize the u c = critical point and will just treat u as a dialable parameter (or just imagine dialing the mass square 
for the vacuum theory). Our main purpose is to identify and understand the critical behavior around the critical 
point which is independent of the specific mechanism to realize it. 

As will be discussed in subsequent sections, as we approach u c = from the uncondensed side (u > 0), the static 
susceptibility remains finite, but develops a cusp at u = and if we naively continue it to u < the susceptibility 
becomes complex. Below we will refer to this critical point as a bifurcating QCP. 



2. Hybridized quantum critical point 

For the second type of instability, one can readily check numerically that is a monotonically increasing function 

of k for negative k+. See fig. 19. Thus the instability goes away if a+\k = 0) > 0, which corresponds to k+ 
satisfying 33 



a 



(o) 



(fc = 0) 



> K + > K C EE ^ . (4.11) 

b { °\k = 0) 
At the critical point k+ = n c 

a ( +\k = 0,K c ) =0 (4.12) 

and as a result the uniform susceptibility \ m (4-5) diverges. Such a quantum critical point was first discussed in [46] 
and elaborated further in [47]. As will be seen below, the presence of the strongly coupled IR sector described by 
AdS 2 gives rise to a variety of new phenomena which cannot be captured by the standard Landau- Ginsburg- Wilson 
paradigm. For reasons to be clear below, we will refer to such a critical point as a hybridized QCP. 



One might be surprised that in the case of a broken continuous symmetry in d = (2 + 1) the ordered phase can persist at any finite 
temperature, as one generally expects fluctuations of the Goldstone mode to destroy the order. As there is only one Goldstone mode (as 
compared to TV 2 other fields) this is a 1/N 2 effect and so is dual to a quantum effect in the bulk; it is shown in [144] that at one-loop 
order in the bulk the symmetry is indeed restored over exponentially long length scales. 

There is also an instability for k+ > 0, but one can check that it corresponds to a UV instability of the vacuum. Here we are interested 
in the IR instability due to finite density effect. See [47] for more discussion. 
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3. A marginal quantum critical point 

We can also tune n + and u together to have a doubly tuned critical point at u = 0, n+ = where the susceptibility 
both diverges and bifurcates. The value of k+ can be obtained from the u — ^ limit of the expression for n c given in 
(4.11), leading to 

< = ~ (4-13) 

where a and f3 are constants defined in Eq. (B17). For the specific example (4.10) of tuning the AdS4 mass of a 
neutral scalar to reach u = 0, the values of a, f3 can be found numerically which give = —2.10. 

As we will see below, the susceptibility around such a critical point coincide with that of the bosonic fluctuations 
underlying the "Marginal Fermi Liquid" postulated in [12] for describing the strange metal region of the high T c 
cuprates. 34 

We note that while understanding the end point of the instability, i.e. the stable phase in which O is condensed, 
requires full nonlinear analysis of bulk equations, to diagnose the instability, to identify the onset of instability (i.e. 
the quantum or finite temperature critical points), and to find the critical behavior near the critical points (from the 
uncondensed side), the linear analysis which gives rise to (4.2) suffices. 

We now proceed to describe the critical behavior near the various critical points identified above. We will use a 
neutral scalar field for illustration. We should also mention that in all critical points identified above, the uncondensed 
side is described by an extremal charged black hole. Given the interpretation in Sec. II F of an extremal charged black 
hole as an intermediate state rather than a genuine ground state, the critical behavior discussed below from the 
uncondensed side should again be interpreted as due to intermediate-energy effects; we return to this point in the 
conclusion. 



C. Quantum Critical behavior: Bifurcating quantum critical point 

1. Critical behavior 

Let us first consider the bifurcating critical point u = 0. Here we will outline the results; for details see the recent 
discussion in [47]. 

To study the critical behavior from the uncondensed side, we take Z/& — » in (4.6), i.e. with both k 2 / ji 2 and u 
small (in particular, below we will always consider k 2 / fi 2 of order or smaller than u which is the interesting regime). 

This limit is a bit subtle as the factor I ^- \ in the SLQL Green function (2.40) behaves differently depending on 

the order in which we take the — » and uo — >• limits. For example, the Taylor expansion of such a term in small 
v\z involves terms of the form log(a;//i*), but in the small uo limit, the resulting large logarithms may invalidate the 
small Vk expansion. We thus proceed carefully. Expanding the explicit form of the AdS2 Green's function (2.41) and 
a±, b± to leading order in z//~, but keeping the full dependence onw, we find that 

k) = xo ) y-Hr{ + ■ ■ ■ (4.14) 

sinh (^log(^)) 



;iv v fi „ . „ ( & 



where 

Xo = nf u ^, u a = 2/i* exp ^ - 7#J , uo b = 2/i* exp ^ - 7^ j , (4.15) 

je is the Euler number, and a,/3,a,/3 are some constants introduced in (B17). 

This expression is exact for all frequencies that are much smaller than the UV scale; however to understand the 
physics, we should understand that the non-commutativity of the uj — ^ and — ^ limits in equation (4.14) 
essentially defines a crossover scale 



Aco ~ exp — , (4.16) 




It has also been pointed out by David Vegh [145] that the retarded function for a scalar operator with v = in AdS2 gives bosonic 
fluctuations of the "Marginal Fermi Liquid" . 
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with # some 0(1) number. 

For uj ^> Kqoi one can expand the arguments of the hyperbolic sine to find 



log fe) 

with the spectral function given by 



log (^-) - if 

xM)=Xo YH + 0(u,k 2 ) (4.17) 



However, for <C Aco, approaching the critical point from the (stable) u > side, we find 



X (uJ 0, fc) = xo - 2i/ fe x* - 4 ^feX* [^~J + • • • • ( 419 ) 

The static susceptibility is obtained from (4.19) by taking — >• 0, and interestingly, near the critical point it 
approaches a finite constant xo fc> r k = 0. It nevertheless develops a branch point singularity at u = 0, as Vk=o — \[^\ 
the susceptibility is trying to bifurcate into the complex plane as we cross u = 0. Upon Fourier transformation to 
coordinate space, this branch point singularity leads to a correlation length that diverges at the critical point. Indeed 
by comparing (4.7) with (2.42), we find that 

(4.20) 



Thus as u — » 0, the correlation length £ diverges as u~i with the same exponent as that in a mean field theory. 

When u < 0, the SLQL scaling dimension of becomes complex for sufficiently small k as = ^Ju + = — 

is now pure imaginary. 35 For a given nonzero Co> and \u\ sufficiently small, the corresponding expression for x(k^oj) 
can be obtained from (4.14) by simply taking = —iXk-, after which we find 

sin (a* log (^)) 

X (u, k) = xo — ) f-Hv + • • • • (4-21) 



sin (A, log 



For uj ^> Aqo one finds the same critical behavior as (4.17), but for uj <C Aco, we find a geometric series of poles in 
the upper-half complex frequency-plane at 

uj n = iuj a exp ^--^L^ ~ iji » n = l,2,... (4.22) 

with 

A /i? = /iexp( ^= j , (4.23) 



indicating that the disordered state is unstable and the scalar operator condenses in the true vacuum. In particular, 
we expect (4.21) to break down for uj ~ uji ~ A/^, the largest among (4.22), and at which scale the physics of 
condensate sets in. This is indeed consistent with an explicit analysis of classical gravity solutions in [45, 47, 139] 
where it was found that O develops an expectation value of order 



Note that the choice of branch of the square root does not matter as (4.14) is an even function of v^. 
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The exponent | in (4.24) is the scaling dimension of O in the SLQL for u = 0, while A is its scaling dimension in 
the vacuum. Furthermore, one in fact finds [45, 47, 139] an infinite tower of excited condensed states in one to one 
correspondence with the poles in (4.22) 

<0>„~M A exp(-^)~M A (^) f n=l,2,.... (4.25) 

This tower of condensed states with geometrically spaced expectation values is reminiscent of Efimov states [146]. 36 
The n = 1 state is the ground state with the lowest free energy which scales as (with that of the disordered state 
being zero) 

F ~ -A IR . (4.26) 




FIG. 20. Left plot: the spirals describing responses of the tower of "Efimov" states. A and B in the figure denote the source and 
expectation value for O respectively. The system under consideration has a Z2 O — »• —O reflection symmetry, which gives rise 
to a pair of infinite spirals related by this symmetry. The condensed states are given by the intersections of the spirals with the 
P-axis with B±^ the ground states and B±^ the first excited states and etc. The red straight line has slope given by (4.27). As 
yj —u — >• most part of the spiral becomes parallel to it. For ease of visulization a nonlinear mapping has been performed along 
the the major and minor axes of the spiral; while the zeros of A and B are preserved by this mapping the location of divergences 
and zeros of -rj are not. Right plot: Response curve in a Landau-Ginsburg mean- field model Clg = \u^ 2 + §(<9^) 2 + A^ 4 + . . . , 
where ip is a Z2 order parameter and the critical point lies at u = 0. The plot is for a small negative u (condensed side). Note 
that the part to the right of behaves as As instead of the linear behavior of the left plot since at the critical point the 
uniform susceptibility diverges. The part between B+ (B-) and P+ (P-) describes a metastable region while the part between 
P+ and P- describes an unstable region. In the left plot the region between P+ and P- is replaced by a pair of infinite spirals 
which give rise to a tower of infinite excited states. 



It is also interesting to compute the response of the system to a static and uniform external source (i.e uo = k = 0) 
in this tower of "Efimov" states, which turns out to be described by a pair of continuous spiral curves given in Fig. 20 
(see [47] for more details). A and B in the figure denote the source and expectation value for O respectively. The 
set of Efimov states are obtained by setting the source A = 0, i.e. they can be read from the intersections of the 
spirals with the £?-axis. For comparison we also show the response curve for the standard Landau-Ginsburg story in 
the right plot of Fig. 20. Approaching the critical point \J—u — >• 0, the spirals exhibit the following features: 

• The spirals are being squeezed into a straight line, with limiting slope 

= Xo (4.27) 
which agrees with the value found from linear response approaching the critical point from the other side (4.19). 



B 
A 



In fact the gravity analysis reduces to a quantum mechanical problem which is exactly the same as that for the Efimov states. 
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• The susceptibility from the condensed side, which is defined by 

XL 



dB 
dA 



however, differs from (4.27), and is given by 

XL = Xo + X* +0(u), 



A=0 



x* = fir 



(4.28) 



(4.29) 



Essentially, even though the spiral is squished into a straight line as we approach the transition, each intersection 
of the spiral with the A = axis has a different slope than the limiting slope of the entire spiral. Note that this 
result is independent of n and in particular applies to n = 1, the ground state. Since xo is the value at u = 
from the uncondensed side, we thus find a jump in the value of uniform susceptibility in crossing u = 0. 

The existence of a tower of "Efimov" states with geometrically spaced expectation values may be considered as 
a consequence of spontaneous breaking of the discrete scaling symmetry of the system. With an imaginary scaling 
exponent, (4.21) exhibits a discrete scaling symmetry with (for k = 0) 

uj^e^^uj (4.30) 

which is, however, broken by the condensate. 37 The tower of "Efimov" states may then be considered as the "Goldstone 
orbit" for this broken discrete symmetry. 

At a finite temperature, using the finite T generalization of (4.6) (as discussed at the end of Sec. II D), equation (4.14) 
generalizes to 



sinh 




log 


(¥) 






) 


sinh 




log 




+*{\ 




) 



X {T) (u,k) = xo- 

where ip is the digamma function. Let us first look at the static regime uj = 0, for which we find 

sinh (v k log (7 



(4.31) 



X {T) (k) = Xo 



sinh (v k log (^)) 



where T a ^ some constants of order fi given by 



T a = e< 

7T 



-t—eP . 

TT 



Similarly to (4.21), the expression for u < is obtained by analytically continuing (4.32) to obtain 

sin (\ k log (|r)) 



X iT) (k)= X o- 



sin 



(A fc log ^ 



(4.32) 



(4.33) 



(4.34) 



Both (4.32) and (4.34) are analytic at u = and for T ^> Aco reduce to 



X {T) (k)= X0 



log^ 
log£ 



0(u,k 2 ) . 



(4.35) 



When T <C Aco, for u > equation (4.32) crosses over to an expression almost identical to (4.19) with uo replaced 
by T. For u < 0, equation (4.34) has poles at (for k = 0) 



T n T a exp 



4/i* 



exp 



a 



-u a 



n e Z+ . 



(4.36) 



37 Note that for n = 1 state, since the physics of the condensate sets in already at Ai, the range of validity for (4.21) is not wide enough 
for the discrete scaling symmetry to be manifest. 
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Comparing to (4.36) and (4.22), we see that T n ~ uo n . The first of these temperature should be interpreted as the 
critical temperature 




(4.37) 



below which the scalar operator condenses and (4.34) is no longer valid. 38 Note that equation (4.34) also has poles 
T n with n < 0. But such T n are of order or much larger than \i for which our approximation is no longer valid, and 
so we ignore them. 

For nonzero cj, in the quantum critical region T ^> Aco, equation (4.31) becomes 




(4.38) 



which can now be applied all the way down to zero frequency. Equation (4.38) reproduces (4.17) for uo ^> T and (4.35) 
for w<T. 

In Fig. 21 we summarize the finite temperature phase diagram. 




FIG. 21. Finite temperature phase diagram with the quantum critical region for bifurcating criticality as a function of u. The 
dotted line is given by Aco in (4.16). But note that the discussions there are not enough to determine the O(l) factor in the 
exponent for Aco- The dynamical susceptibility in the bowl-shaped quantum critical region is given by (4.38) with the zero 
temperature limit given by (4.17). 



2. Physical interpretation: SLQL as a fractionalized phase of Bose condensates 

Now let us now try to interpret the above results. First we emphasize that nowhere on the uncondensed side 
do we see a coherent and gapless quasiparticle pole in the dynamical susceptibility, which usually appears close to 
a quantum phase transition and indicates the presence of soft order parameter fluctuations. The facts that at the 
critical point the susceptibility (4.19) does not diverge and the spectral function (4.18) is logarithmically suppressed 
at small frequencies are also manifestations of the lack of soft order parameter fluctuations. 

It has been argued in [44, 47] that the quantum phase transition at a bifurcating critical point corresponds to 
a confinement /deconfinement transition. Across the critical point, the scalar operator Ok=o develops a complex 



Including frequency dependence, one can check that as T is decreased through each T n , a pole of x^ T )(o;, k) moves from the lower half 
o;-plane to the upper half-plane. Thus we see the interpretation of each of these T n ; as the temperature is decreased through each of 
them, one more pole moves through to the upper half-plane. There exist an infinite number of such temperatures with an accumulation 
point at T = 0; and indeed at strictly zero temperature there is an infinite number of poles in the upper half-plane, as seen earlier 
in (4.22). Of course in practice once the first pole moves through to the upper half-plane at T c =T±, the uncondensed phase is unstable 
and we should study the system in its condensed phase. One can further study the critical behavior near the finite temperature critical 
point T c . Here one finds mean field behavior. 
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dimension in the SLQL and conformality is lost. The loss of conformality is realized through a dynamically generated 
"confinement" scale Air below which an infinite tower of geometrically separated bound states of operator O form and 
then Bose condense. The physical picture here is similar to the BEC regime in a strongly interacting ultracold Fermi 
system where fermions form bound molecules and then Bose condense. The bifurcating QCP is thus characterized by 
the onset of forming bound states rather than by soft order parameter fluctuations; this may help explain why the 
susceptibility does not diverge at the critical point and the spectral density is suppressed. This should be contrasted 
with the fermionic case discussed in Sec. HIE, where fermionic bound states each form a Fermi surface. In both cases, 
SLQL serves as a deconfined high energy state for the more stable lower energy states which arise from SLQL via 
bound state formation. 

To summarize, for a bifurcating QCP, while the phase transition can still be characterized by an order parameter, 
the order parameter remains gapped at the critical point and the phase transition is not driven by its fluctuations. 
Instead the phase transition appears to be driven by formation of bound states. 



D. Quantum critical behavior: a hybridized QCP 



1. From uncondensed side 



We now examine the critical behavior around a hybridized QCP (4.12) (again using a neutral scalar) following [46, 
47, 142]. 

The full dynamical behavior from the uncondensed side can be obtained by expanding a+(cj, fc) in (4.6) around 
uj = 0, k = and k+ = k c , which leads to 



1 



n+ — k c + h k k 2 — h^uj 2 + hQk(oj) 



with 



hh = 



«<?>(*) 

it'll) — 



fe=0,, 



h = 



2 "*o<_ 0) (fc) 



k=0 



K + = 



6f)(Jk) 



(4.39) 



(4.40) 



k=0. 



Recall that the SLQL retarded function Gk(w) ~ uo 2vk . From explicit gravity calculation one finds that the various 
constants in (4.39) have the following behavior: h k > 0, h < 0, and h u > (for Vk=o > !)• 
Setting uj = in (4.39), the static susceptibility x(fc) can be written as 



X{k) « — 

(«+ - k c ) + h k k 2 

which has a standard mean field behavior, with the spatial correlation length scaling as 



f~ (K+-K, c ) "crit = ^ ' 

and critical exponents (see Appendix F for a review of definition of critical exponents) 

7 = 1, r] = . 



(4.41) 



(4.42) 



(4.43) 



The behavior of the full dynamical susceptibility (4.39) depends on the competition between the analytic contribu- 
tion h^u 2 and the non-analytic contribution Gk{w) from SLQL. When v^ = q G (0, 1), the non-analytic part dominates 
at low energies and the analytic contribution can be ignored, leading to the interesting and decidedly non-mean-field 
behavior 



(«+ - k c ) + h k k 2 + hC(v)(-iuj) 2u 



(4.44) 



where v = v k =§ and constant C(y) is real and negative for v G (0, 1). We will consider k = below. At the critical 
point k+ = k c we find that 



X(u,k = 0) - (~ioj) 



-2v 



(4.45) 
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Away from the critical point, the relative magnitude of the two terms (with k = 0) in the denominator of (4.44) 
defines a crossover energy scale 

A&>~|« c -/t + |& . (4.46) 

For u < A ( co one finds (4.41) while for u > (4.45). This crossover scale (4.46) defines the correlation time £ r 
of the system 

A co 

Comparing the above expression with (4.42) we then find that £ r ~ £ z with the dynamical exponent z given by 

z = i . (4.48) 

Of course this exponent can equivalently be seen by balancing the k 2 term and the uj 2v term in (4.44). Also note that 
when n+ < n c equation (4.44) has a pole in the upper half plane (since hC(y) > 0), 

^pole ~ iA ( cO ( 4 « 49 ) 

indicating an instability. 

When v > 1, in (4.39), the non-analytic part Qk(w) ~ u 2u from SLQL becomes less important than the analytic 
corrections ~ u 2 . Since the term proportional to Qk(w) is the only complex term, it should nevertheless be kept. 
Now (4.39) describes a long-lived (nearly gapless) relativistic particle with a small width T ~ uj 2v '. The dynamical 
exponent is now z = 1. 

At a finite temperature T <C /i, equation (4.39) generalizes at leading order in T/fi to (see discussion at the end of 
Sec. II D), 

X(u, T) » ^ =y— (4.50) 

where /i^T a constant) 39 comes from (analytic) finite temperature corrections to a + and b+. The finite tem- 
perature SLQL retarded function Q^\uj) has the form Q^\uo) = T 2vh gb(^,Vk) with a universal scaling func- 
tion (see (2.49) and (2.50)). Consider first uj = in (4.50), which is 

v( T ) 1 . (4.51) 

It is interesting that the analytic contribution now dominates for v > \. For v < | we find that there is a pole at 

T c -(^ c -^+)A (4.52) 

for k + < k c . This should be interpreted as the critical temperature for a thermal phase transition, above which the 
instability disappears. From the uncondensed side, such a temperature scale coincides with the crossover scale (4.46). 
For v > |, we find instead mean field behavior 

T c ~ (k c - «+) (4.53) 
and the finite temperature crossover scale becomes 

A^~|« c -« + | (4.54) 

which no longer tracks the scale (4.46) of zero temperature. With nonzero lj, when v < |, we find uj/T scaling, while 
for v > \ there is no uj/T scaling. 



From explicit gravity calculation one finds Ht > for iSk=o > i. 
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2. Low energy effective action 

The discussion in Sec. HE 2 leading to the low energy effective action (2.75) for the fermionic case can be applied 
almost identically to a scalar field. In particular, near a hybridized QCP (4.12), K in (2.72) can be expanded near 
uo = and k = as (say for a neutral scalar) 

K= a+pk) = +fc p_^ w 2 + _. (455) 

b+(u,k) 

We thus find that near a hybridized QCP, the low energy effective theory can be written as 

S eJ f [ip,<f>} = Sslql [*] + A j 'dt&/> + S LG + j 'ipj (4.56) 
where Sslql is the same as (2.77) (except that is now a scalar operator) and 

Slg = -\J ^K(k,uj)xp (4.57) 

can be considered as the quadratic part of the Landau-Ginsburg action for a scalar order parameter ip near criticality. 
The full low energy action (4.56) describes the hybridization between a Landau-Ginsburg sector and a strongly coupled 
sector controlled by SLQL. This mixing between two sectors is the key to understanding the critical physics, and is 
why we call this a hybridized critical point. 

The behavior of the dynamical susceptibility at T = discussed earlier and in particular the crossover to mean 
field behavior at v = 1 can be readily seen from the effective action (4.56). Recall that <$> has SLQL scaling dimension 
| + z/, and from the second term in (4.56), ip has dimension - — v. For v > 1, the dimension for ip become smaller 
than — |, and as a result the kinetic term (dtip) 2 (coming from uJ 1 term in (4.55)) in (4.57) becomes relevant and 
more important than the hybridization term &ip (which is marginal by definition) in (4.56). Alternatively, we can 
now assign — | as dimension of ip using (dtip) 2 , under which the hybridization term $>ip will have dimension v which 
becomes irrelevant for v > 1. 

The crossover to mean field behavior at v = \ at a finite temperature discussed in (4.51)-(4.53) can also be readily 
seen from the effective action (4.56). Finite temperature can be thought to generate a term J dtTip 2 , which becomes 
relevant when the dimension of ip becomes smaller than zero, i.e. for v > |. Alternatively we can now use Tip 2 term 
to assign dimension to ip, under which the hybridization term <&ip then becomes irrelevant for v > ^ as now the 
dimension for <l> becomes larger than 1. 

Finally, we mention that for v < ^ there is an alternative description in which the SLQL theory is studied in its 
alternative quantization, i.e. a different IR operator <£_ is used that has dimension \ — v. We do not go into details of 
this description here except to point out that in this description the critical point ft + = k c can be understood as tuning 
the relevant deformation <& 2 _ in the SLQL to 0. This is thus consistent with the standard Landau-Ginsburg- Wilson 
paradigm of a critical point as a quantum field theory with a single relevant direction; however the interaction of this 
relevant direction with the sector described above by the field ip is nontrivial, resulting in the physics described above. 
We refer the reader to [46, 47] for more discussion of this. 



3. Free energy: condensed side 



Finally let us look at the the scaling of the expectation value and the free energy from the condensed side. This 
can be done by analyzing the condensed solution on the gravity side [46] . Alternatively, one could use the low energy 
effective action (4.56) [47, 142] which provides a significant simplification. We will use the latter approach. We first 
generalize (4.56) by including the next nonlinear term for ip in (4.57), i.e. 

S lg = ~\ j ^K(k,uj)iP-u J dtip 4 + ... (4.58) 

with u some positive constant. 

Now consider that ip develops some nonzero expectation value. From the relative scaling dimensions between <1> 
and i/j, we can then write & as 



(4.59) 
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and the free energy density F corresponding to (4.56) can then be written as 



F ~ Cip + -(«+ - k c )^) 2 + ^ 4 



(4.60) 



where the first term comes from the Qip term with C some constant. Equation (4.60) can also be derived from a 
detailed bulk analysis which also gives that C > for v < \. Now notice that for ip small, the first term dominates 
over the ip 4 term if v < |, while the Landau-Ginsburg ip 4 term dominates for v > |. In other words, since the first 
term is marginal by assignment, ip 4 term becomes relevant when v > ^. 40 
For v < | we can ignore the last term in (4.60) and for ft + < k c find that 



and as a result 



(4.61) 
(4.62) 

(4.63) 

where in the second equality we have expressed the free energy in terms of correlation length using (4.42). From (4.63) 
we can also deduce that at the critical point we should have 



Including the source J, which has dimension \ + ^, the free energy should then be given by a scaling function 



(O) ~iP~ Jh + " . 



(4.64) 



From (4.61), (4.64) and (4.63) we can collect the values of various scaling exponents (see Appendix F for a review of 
exponents) 



« = 2- * 0=ip, 5=^ 
2u 2v k — v 



(4.65) 



For v > |, we can ignore the first term in (4.60) and the analysis becomes the standard Landau-Ginsburg one. As 
a result, the behavior near the critical point becomes those of mean field [46]. We thus find that for v > |, 



(0) ~ ^ ~ (k c - F ~ -(« c - k + ) 2 (0) K+=Ke ~ J* 

and various exponents become 



a = 0, /?=-, £ = 3 



(4.66) 
(4.67) 



which agree with the values of (4.65) for v = \. 



4- Summary 

We summarize the somewhat intricate critical behavior near a hybridized QCP for various values of v in the following 
table 



Quantity 


" e (o, \) 




^a,i) 


v > 1 


uo/T scaling 


yes 


yes 


no 


no 


A M z 


HYB 


HYB 


HYB 


MFT 


T c 


HYB 


HYB 


MFT 


MFT 


a, /3, 5 


HYB 


MFT 


MFT 


MFT 


1 1 Vi ^crit 


MFT 


MFT 


MFT 


MFT 



40 Some readers might worry that higher powers like ip 6 may also become relevant at some point (for example for v > |). But note that 
once the last two terms in (4.60) dominate we should reassign the dimension of ip and the standard Landau-Ginsburg story applies. 
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In the above "MFT" denotes mean-field behavior and "HYB" (which is for hybridized) refers to the following non- 
mean-field exponents for zero-T crossover scale A^q, the dynamical exponent z and critical temperature T c , 

A&£~(« c -«+)^, s=i T c ^(^ c -^ + )^ (4.68) 



and 

1 _ l 

2 

2v 1 2^ 



For all values of v the static susceptibility from the uncondensed side is always given by the mean field behavior (4.41). 
Note in particular that for v € (0, |) one finds an expression 

X(u, k; T) » i (4.70) 

K+-K c + /l fc fc 2 +/lT2^ b (^) 

with ^5 given by (2.50) (with q* = 0). Equation (4.70) is of the form (1.17) for the dynamical susceptibility of 

CeCuQ- x Au x with a choice of 2v = 0.75 [46], although the universal function appears to be different from the 
f(u/T) there obtained by fitting experimental data. 



E. Doubly fine-tuning to a marginal critical point 

We have just described two different kinds of fixed point; the bifurcating transition involves tuning u — » 0, and 
the hybridized transition involves tuning n+ — >• k c . We can obtain a new kind of critical point by tuning u — >• and 
n+ — » k c at the same time, at some critical value k* given by 

< = ~ (4-71) 

where we have used the v limit of the expression for k c given in (4.11) and a, /3 are defined in (B15)-(B17). We 
will fix u — and vary k+. 

The dynamical susceptibility can be obtained from the ac + = result at the bifurcating critical point (4.17) and (C6) 



log (^) -i% 

X^H^k) = If 7 / - f s V . (4.72) 



Now using the definitions of uo a ^ in (4.15), applying the Wronskian relation (B20), and expanding near ~ ft* = — ^, 
we find: 



log(^) - if 

X ( * +) (",*) = ? V% \ — • (4.73) 



(«+ - ft c *) (log (- 



This contains a good deal of interesting physics. First let us keep k+ > and study the very low-frequency behavior: 
expanding in powers of the inverse logarithm, we find 

x («+)( w ->0 ) fc)~ ( 4 - 74 ) 

K+-K* \logUjJ 

Thus the static susceptibility diverges as k —> «*, just as for the hybridized critical point discussed above. Away from 
the critical point the leading contribution to the spectral density appears at order log - (w), 

Imx (K+) (w^0,fc)~- . „ log" 2 ( — ) +o(—^-) (4.75) 

4(/t+-/t;K W Viog 3 ^/ 

This is the same as the corresponding result close to a bifurcating critical point with zero double-trace deformation. 
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Next, we note that (4.73) contains a pole in the upper half plane at 



cj* = iu b exp ( 1 1 — - j (4.76) 

Note however that if ft + > ft*, the pole is at an exponentially high energy scale, this low- frequency formula breaks 
down far before the pole, and the pole is not physical. This is the stable phase. However, if ft + < ft*, then close to the 
critical point the pole is very close to the origin and this formula is valid. The instability pole is now physical, and 
we see that a new low-energy scale has been exponentially generated. This is reminiscent of the effect of a marginal 
coupling growing strong in the infrared, as in QCD or the BCS instability. Indeed as discussed in (2.59) and (2.60), we 
expect in the IR a double trace deformation to SLQL will be generated. At u = = k : the double trace deformation 

5S= ' 



\1 -TL) n 

*%-Jdt&, (4.77) 



has dimension 1 and is thus marginal It can be checked explicitly that the value of the effective marginal double-trace 
parameter n\n changes sign at ft + = ft* (although one must be somewhat careful; for details see [47]). Thus what 
we have found above is that if ft + > ft* then the coupling is marginally irrelevant and has no effect in the IR, whereas 
if ft + < ft* it is marginally relevant and drives the theory to a new fixed point (i.e. the scalar will condense). 

Similarly, one can compute the corresponding finite temperature correlator using (4.38); a precisely analogous 
calculation gives 

X (K+ \u,k;T) = , \ b \ • (4-78) 

(«+ " «S) ( lo S (^r) i£r)) ~ 27ZW 

From here we can see that for ft + < ft* the static susceptibility diverges at the critical temperature 

T c = T b exp ( — (4.79) 

\2vuj3 2 ft + - ftjy 

above which the system is stable. The temperature is set by the same dynamically generated scale. 

Now we consider the fluctuation spectrum at criticality. Returning to T = 0, we now keep uj ^ and sit precisely 
at the critical point ft = ft*. From (4.73) we now find 

X (k+=k:) (u, k) = -2v v p 2 {log (j^j - i^j (4.80) 

The appearance of a pure logarithm in the numerator of this propagator at criticality is interesting. Relatedly, at 
criticality the spectral density 

Im X («+=«:) = 7r^/3 2 sgn(cj) (4.81) 

is a pure step function 41 . This should be contrasted with the situation away from criticality, (4.75), in which there 

is a logarithmic suppression at low frequencies; as expected, by sitting at criticality we find significantly more low 
frequency modes. 

Similarly, setting ft — >• ft* in the finite temperature expression (4.78) we find 

X^ += <\uj,k;T) = -2v v (3 2 (log Q "^)) ( 4 - 82 ) 
Taking the imaginary part we find 

Im X («+=<)( w ,A ; ;r)= 7 r I /a/3 2 tanh(^) (4.83) 



The logarithm jumps by in as we cross through uo = 0, resulting in the step function; note that this was necessary in order to maintain 
the relation o;Im x{ u ) > 0> true for any bosonic spectral density. 
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which is simply a smoothed-out version of the step function that we find at zero temperature. See Fig. 22 for the 
finite temperature phase diagram for a marginal critical point. 

This form of the finite-temperature fluctuation spectrum is potentially of great physical interest. Note that (4.83) 
implies that 

Tm X (« + =<)(u,,k;T)~\% "^ T ^ T (4.84) 
[sgn(cj) u > T 

which is precisely of the form for spin and charge fluctuations in the phenomenological "Marginal Fermi liquid" [12] 
description of High-T c cuprates in the strange metal region (as discussed around equation (1.18) in Sec. IC). Thus 
this marginal critical point can be viewed as a concrete realization of the bosonic fluctuation spectrum needed to 
support a Marginal Fermi liquid. 




FIG. 22. Finite temperature phase diagram with the quantum critical region for marginal criticality at u — and changing 
— «+). The susceptibility in the bowl-shaped quantum critical region is given by (4.82) with the uj ^> T limit given 
by (4.80) 



F. Summary 

In this section we have discussed several types of quantum critical points related to SLQL which to different degrees 
lie outside the Landau- Ginsburg- Wilson paradigm: 

• A hybridized QCP describes an order parameter if; with a Landau-Ginsburg effective action Slg hybridized 
with an operator <1> in SLQL. Depending on the scaling dimension of <1> in the SLQL, the phase transition could 
exhibit a rich spectrum of critical behavior, including locally quantum critical behavior with nontrivial uj/T 
scaling. At the level of effective theory, this critical point lies mildly outside the standard Landau paradigm, 
as the phase transition is still driven by soft fluctuations of the order parameter and all the critical behavior is 
fully captured by (4.56). As discussed in Sec. HE, on the gravity side the Landau-Ginsburg sector is associated 
with the appearance of certain scalar hair in the black hole geometry, which lies outside the AdS2 region. 

• A bifurcating QCP arises from instabilities of the SLQL itself to a confined state and is not driven by soft order 
parameter fluctuations. On the condensed side, a scalar operator develops a complex scaling dimension in the 
SLQL, generating a tower of bound states, which then Bose-Einstein condense (the possible bound-states form 
a geometric series of exponentially generated scales). 42 In particular, one finds a finite critical susceptibility 
with a branch point singularity, and the response of condensed states is described by an infinite spiral. As 



SLQL may be considered as a "deconfmed" state in which the composite bound states deconfme and fractionalize into more fundamental 
degrees of freedom. 
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discussed in [47], at a field theoretical level, underlying these features is the annihilation (and moving to the 
complex plane of a coupling constant) of two conformal fixed points, which is very different from the standard 
Landau- Ginsburg- Wilson paradigm of phase transitions, which can be characterized as a single critical CFT 
with some relevant directions. On the gravity side, for charged operators the instabilities of the SLQL can be 
understood as the pair production of charged particles which then subsequently backreact on the geometry. For 
a neutral scalar operator, the instability is related to the violation of the BF bound in the AdS2 region. 

• A marginal QCP can be obtained by sitting at the critical point of a bifurcating QCP and then dialing the 
external parameter which drives a hybridized QCP. Here given the critical theory describing a bifurcating QCP 
(which comes from the merger of two fixed points), the phase transition can be described as the appearance of 
a marginally relevant operator. Interestingly, the fluctuation spectrum that emerges is identical to the bosonic 
fluctuation spectrum that (when coupled to a Fermi surface) is postulated to underly the "Marginal Fermi 
Liquid" description of the optimally doped cuprates [12], making this critical point of potential importance. 
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FIG. 23. How to interpret the quantum critical behavior of this section. SLQL should be interpreted as a universal intermediate 
phase which orders into some other phases, such as a Fermi liquid, at lower energies. The dome in the figure outlines a region 
near the critical point to which our current discussion does not have access. The discussion of this section only describes the 
quantum critical behavior outside the dome-shaped region. 

Finally let us elaborate on an important point, which we have glossed over in our discussion so far. As discussed in 
Sec. II F, SLQL, which describes the disordered phase in our examples above, should be interpreted an intermediate- 
energy phase, rather than a genuine ground state. That is, we expect SLQL to order into some other phases at lower 
energies, which may not be visible at the large N limit we are working with. An example discussed in Sec. HIE is 
that at some exponentially small scale in TV 2 , SLQL orders into a Fermi liquid phase. 43 Thus the quantum critical 
behavior found in this paper should be more correctly interpreted as describing the intermediate-region indicated in 
Fig. 23. 



V. DISCUSSION AND CONCLUSIONS 



In these lectures, we used gauge/gravity duality to study a class of finite density systems, which are described on 
the gravity side by a charged black hole. We showed that the duality predicts a universal intermediate-energy phase, 
called a semi-local quantum liquid (SLQL). Such an unstable phase is characterized by a finite spatial correlation 
length, but an infinite correlation time and associated nontrivial scaling behavior in the time direction, as well as a 
nonzero entropy density. We discussed two sets of phenomena related to the SLQL: 



See [143] for a recent discussion of nucleation of a neutral order parameter in a Fermi liquid-like phase. 
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1. Gapless fermionic or bosonic degrees of freedom (described by mean field) hybridized with degrees of freedom 
in SLQL. In either case, the system can be described by the following low energy effective theory 



In the fermionic case, ^ denotes excitations around a Fermi surface with S'mean field [^] given by that for a 
free fermion (plus possibly higher order self-interactions). In the bosonic case, \I/ denotes the order parameter 
around a quantum critical point with S mean field [^] its Landau- Ginsburg effective theory. In both cases, is some 
(fermionic or bosonic) operator whose dynamics is controlled by the strongly coupled SLQL. The IR behavior of 
the theory (i.e. whether one finds a Fermi surface with or without quasiparticles, or a quantum phase transition 
which obeys locally quantum critical or mean-field behavior) depends crucially on the scaling dimension of in 
the SLQL. With a suitable choice of the SLQL scaling dimension, IR phenomena which strongly resemble those 
of high-T c cuprates in the strange metal phase (including a linear resistivity) and that for CeCuQ- x Au x near a 
quantum critical point were found. 

On the gravity side the gapless degrees of freedom are associated with the appearance of certain scalar or 
fermionic hair in the black hole geometry, and lie outside the AdS2 region which describes the SLQL. 

2. Instabilities of the SLQL itself. We discussed examples in which SLQL orders at lower energies to a Fermi 
liquid, superconductor, or antiferromagnetic-type state. In all cases certain operators develop a complex scaling 
dimension in the SLQL generating a tower of bound states. In the fermion case, each of these bound states forms 
a Fermi surface. In the scalar case, they form a tower of Bose-Einstein condensed states. In particular, near 
the bifurcating critical point where the scalar instability sets in one finds a finite susceptibility and the response 
is described by an infinite spiral. In all of these examples the lower energy states have no zero-temperature 
entropy density and SLQL may be considered as a "deconfined" state for them in which the composite bound 
states deconfine and fractionalize into more fundamental degrees of freedom. This may help explain the nonzero 
entropy density for the SLQL. As discussed in Sec. II F, we expect this picture to be generic; while the precise 
nature of the lower energy state depends on the specific dynamics of the individual system, the SLQL emerges 
universally at intermediate energies through fractionalization. 

Instabilities of SLQL in turn provide novel mechanisms for the emergence of Fermi liquids, superconductors and 
AFM-type of states at low energies. 

On the gravity side, for charged operators the instabilities of the SLQL manifest themselves as pair production 
of charged particles which then subsequently backreact on the geometry. For a neutral scalar operator, the 
instability is related to the violation of the BF bound in the AdS2 region. 

It is natural to ask whether some phase similar to the SLQL could underlie some known strongly correlated 
condensed matter systems, and in particular, whether the observed non-Fermi liquid behavior and novel supercon- 
ductivities in various systems could be attributed to a similar intermediate-energy phase. 

When one encounters scaling behavior in an observable, an important immediate question is whether the behavior 
is due to intermediate-energy or vacuum effects. Here we provide an explicit example in which Marginal Fermi Liquid 
behavior - such as that proposed for high-T c cuprates in the strange metal phase (including linear resistivity) - 
can arise from intermediate energy effects. Local quantum critical behavior similar to that for CeCuQ- x Au x near a 
quantum critical point also appears here as due to intermediate energy effects. 

In many heavy electron systems, while quantum fluctuations from the quantum critical point corresponding to the 
onset of magnetism provide a natural starting point for understanding the observed non-Fermi liquid behavior and 
sometimes novel superconductivity, however, such exotic behavior does not always arise in the proximity of a quantum 
critical point (see e.g. [10]). There thus appears plenty of room for them to arise from an intermediate-energy phase. 

More generally, we expect that candidates for the SLQL to occur include systems which exhibit frustrated or 
competing interaction terms in their Hamiltonian. Such systems can have a large number of near-degenerate states 
near the vacuum, similarly to the holographic systems considered here. Also systems which involve strong competition 
between tendencies towards itinerancy and localization could exhibit the semi- local behavior found here. As mentioned 
earlier, there are also some tantalizing parallels between properties of the SLQL phase and models from dynamical 
mean field theory (DMFT) [33], such as scaling in the time direction with spatial directions as spectators, and a finite 
entropy. The fact that DMFT techniques have been very successful in treating many materials in some intermediate 
energy region provides some comfort in thinking that an underlying universal intermediate-energy phase like the SLQL 
could be at work. 
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Appendix A: Operator dimensions and retarded functions in AdS2 

Here we give a derivation of the scaling dimension (2.37)-(2.38) and retarded Green's function Gk(w) (2.40) for a 
charged scalar in AdS2 x The background metric and gauge field is given by (2.30) which we copy here for 

convenience 



ds 2 = -£(-dt 2 + d( 2 ) + £R 2 dx 2 , A t = -y. (Al) 

We consider the following quadratic scalar action 

S = ~ J d d+1 x^g~ [g MN (d M + iqA M )<P*(d N - iqA N )<t> + m 2 <j>*<f\ . (A2) 

Expanding <fi in terms of Fourier modes in spatial directions, 

dk 



4>(t,x,() = J 7^u^^(t,C) (A3) 
equation (A2) reduces to an action for <j)^ in AdS2 

S = - j d 2 x^-g [g ab (d a + iqA a )<f>t(d b - iqA b )4> R + mj^g] (A4) 



with a mass square 



m^^-^+m 2 , k 2 = \k 2 \. (A5) 

niR 



The indices a, b only run over t, (. Further going to frequency space one finds that equation of motion from (A4) can 
be written as (for notational simplicity we will just denote 4> k (uo, Q as (j) below) 

-d 2 c + V(()0 = O (A6) 

with 

2d2 / „ \ 2 Q/ _ 9 _ 1 



m k R 2 ( , Q*\ _ 2 2uj q* v k 



V(C) = ^-^+^J = _^_^ + ^__a (A7 ) 



where 



1 k 2 R 2 
4 ~ V"" ^ ' 4 + 



v k = \\m\R 2 -q 2 + - = \\m 2 R 2 - ql + - + -^f , q* = qe d . (A8) 



Approaching the AdS2 boundary ( — >> 0, equation (A6) has general solution 

<t> % {u, C -> 0) - a(i~ Uk + P(^ Uk . (A9) 

Compare (A9) to (2.7); it implies that the conformal dimension of the operator to ^ has dimension 5 k = \ + Vk-> 
as claimed in (2.37) and (2.38). 

Introducing x = —2iuj(, equation (A6) can be written as 

^+(4+f +i ^V= o ( ai °) 
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which is precisely the Whittaker's equation and the two linearly independent solutions can be written as 

(j> = C in W iq ^ Uk (x) + C out W- iq ^ Uk (-X) , (All) 

where W\ itl (z) is the Whittaker function. Of these, the function multiplying ci n , 

W Wfc (-2iuj() = e^(-2iujQ^ + . . . , ( -+ +00 (A12) 
is ingoing at the horizon. Keeping only the ingoing part, the asymptotic behavior of <j) at £ — » is given by 

- - r(2 ^ } -(-2^0*-^ + /I r( " 2l/fc) . A -2iu,Ql+»» . (A13) 



r(| + z/fc-ig*) r(|-z/ fe -%) 

From (2.11) the retarded function is simply the ratio of two coefficients in the above equation, leading to 



which is (2.40). 



T(2v k )T \2~ v k- iq*) 



Appendix B: Derivation of the master formula 



Here we present the derivation of the master formula (2.56), following [40]. The equation of motion (in momentum 
space) for with action (A2) in the charged black hole geometry (2.14) can be written as 



zd+ld * [J^ 9 ^) + z " ( (UJ + q f At) ~k 2 j<P- m 2 R 2 <p = . (Bl) 

We will consider T — with / and A t given by (2.26). Recall from (2.7) that near the AdS^+i boundary <j)(z) has 
the standard asymptotic behavior <p(z — ^ 0) ~ Az d ~ A+ + Bz A+ and that the correlator we seek to compute can be 
written as 44 

G R (w,k) = j (B2) 

provided that <j) is an in-falling wave at the horizon. Note that since A t — >• \i as z — >• 0, the boundary theory energy 
corresponds to uj + g/i, i.e. uj should be interpreted as the measured from the effective chemical potential q\i (which 
we take to be positive for definiteness) . 

We are interested in the behavior of Gr in the low frequency limit {uj <C /i). However, we cannot directly perform 
a perturbative expansion in uj in (Bl) since the cj-dependent term becomes singular and dominates at the horizon 
where f(z*) — >• 0, no matter how small uj is. This is not just a technicality; physically this reflects the abundance of 
critical modes coming from the IR fixed point and is of tremendous importance. 

To deal with it, it is useful to isolate the near- horizon AdS2 x region (2.30) (which will be referred to as the 

IR region) from the rest of the black hole spacetime (which will be referred to as the UV region). In the IR AdS2 
region the uj dependence will be treated non-perturbatively. The crossover between the two regions happens near the 
boundary of AdS2 with ^ <C C <^ ^? where fi( ~ remains big so that we are still in the near-horizon region, but 

the cj-dependent term z*-jr ~ UJ 2 ( 2 in (Bl) becomes small. In the UV region we can treat the small uj limit using a 
standard perturbative expansion with the lowest order equation given by setting uj = in (Bl). 

So we start with a solution in the AdS2 x region and evolve it outwards, matching it to a solution in the UV 

region to determine the coefficients A and B. The scalar wave equation on the AdS2 x ^ d ~ 1 region (2.30) is given 
by (A6). As we are computing a retarded correlator we are looking for an in-falling solution at the black hole horizon 
(i.e. the term proportional to Q n in (All) even though the explicit form is not important for the current discussion). 
We now take this infalling solution and evolve it to the boundary of the AdS2 region, where it can be written as 

m=C i-^ + g k{oj)C h+^. (B3) 



We will consider the standard quantization in our discussion below. The generalization to the alternative quantization is immediate. 
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In (B3) we have normalized the solution so that the first term has coefficient 1 and then by definition the coefficient 
of the second term is precisely the retarded Green function in the AdS2 region we computed earlier in (A14). 

We now need to continue evolution all the way to the asymptotically AdS^+i region to extract the coefficients A 
and B. We thus need to match the IR solution (B3) to a solution in the UV region. Now that we are out of the 
dangerous IR region, to leading order we can set uj = in (Bl) in the UV region, 



d z 



f 



yd-l 



d z 



f 



m 2 R 2 d> = . 



(B4) 



The above equation has two independent solutions rj±^ which can be specified by their behavior near z 



z^ as 



»/£ 0) (*)->! C 



with the corresponding asymptotic behavior as z — >• as 



(z — » z*) 



(0)/ 



d-A 4 



(B5) 



(B6) 



a±\k) and b±\k) thus defined are (dimensionless) functions of k which can be computed numerically. 
In the overlapping region both (B3) and (B5) apply, which determines the full UV solution to be 

cj ) (z)=r ] ( °\z)+g k (oj)z 2 ^r 1 ^\z) . 
Equation (B7) can be generalized to higher orders in uj (for details see [40]) 

where 

V±=V ( ± +"V± +0(u 2 ) 



(B7) 



(B8) 



(B9) 



are the two linearly independent perturbative solutions to the full UV region equation. Near z = 0, rj± have the 
expansion of the form (B6) with various coefficients a±\b±* replaced by a±,b± which have an analytic cj-expansion 
such as 



a + (k,uj) =a { +\k)+ua { l\k) + ... . 



(BIO) 



Note that since both the boundary conditions specifying rj± (B5) and the equation (Bl) are real, a±,b± are real. 
From (B8) and the expansion of r]± near z = we thus find the boundary theory Green's function to be 



a + (o;, k) + a_(cj, k)<3k(u)^* k 



(Bll) 



We conclude this discussion with some remarks: 



1. At uj = 0, the leading order equation (B4) in the outer region is in fact the full equation of motion; there is no 

IR region. rj±^ of (B5) now provide a set of basis for the full equation of motion. Note that by construction 77^ 

is normalizable at the horizon. When vanishes, from (B6) rj^ becomes normalizable also at the boundary 
and thus is now a genuine normalizable mode in the black hole geometry. 

2. Note that for a neutral scalar with q = 0, equation (Bl) only depends on uj 2 and the expansion parameter 
in (B9) and (B10) should be uj 2 , i.e. 

a + (k,uj) = af(k)+u 2 a { +\k) + 0(uj a ) (B12) 



and so on. 
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3. The functions a±(cj,k), b±(u,k) are obtained by solving equation (Bl) perturbatively in uj in the UV region. 
Their ^-dependence comes from two sources, from dependence on via the boundary condition (B5) and from 
k 2 dependence in the equation (Bl) itself. Since the geometry is smooth throughout the UV region we expect 
the dependence on both and k 2 to be analytic. In fact we can think of b± and a± as functions of i.e. 
there exists a function b(i/k, & 2 ,o;), analytic in all its arguments, such that b± = b(±Vk : k 2 ,uj). This is clear from 
the boundary condition (B5) (and its generalization for higher orders in uj) and from the fact that there is no 
other dependence on from the equation of motion itself. 

4. As Vk 0, the two solutions in (B5) become degenerate, i.e. 

a+(fc, uj) a_(fc, lj), b+(k : uj) b-(k : uj) . (B13) 

At Vk = 0, the basis of functions in (B5) should be replaced by 

i i 
v (0) {z) = ( 1) \ fj(0) {z) = -(±-Y log ^ (B14) 

where the asymptotic behavior for them at z —> is 



,<., (2) , a (£) S - 4 + ,(£) 4 , *»M-a(±)~ + *(±) 4 (B15) 

o;, /3, a, f3 are now dimensionless real numbers which can again be found numerically. More explicitly, in the 
limit of v\z — >• the basis of functions (B5) can be expanded as 

n f = rf°\z) ± W 0) M + 0(vl) (B16) 

which leads to 

bf = ± v k j3 + . . . , <4 0) = a db i/ fc a + . . . . (B17) 

5. Coefficients a±,6± are not independent. For example evaluating the Wronskian of (Bl) (for uj = 0) 45 for rj±^ 
and demanding that it be equal at infinity and at the horizon, we find the elegant relation: 

(k)b {0) (k) - a (0) (k)bf (k) = — . (B19) 



A similar analysis on 77, fj results in 



af3-f3a = . (B20) 



6. When = —i\k is pure imaginary (i.e. when (2.79) is satisfied), and the basis of solutions (B5) now has the 
form 

These boundary conditions are now complex, and thus so are the rj±. As the rj± actually obey a real wave 
equation, the full analytic structure is determined by the boundary conditions in the infrared; thus we find that 
now t? + = 77* . This also implies that a± , b± are complex and 

a*+ = a_, b\=b_ . (B22) 

7. The derivation leading to (Bll) can be immediately generalized to a finite temperature T <C fi. In the IR 
region one replaces the equation (A6) by the corresponding equation for the AdS2 black hole (2.33). As a result, 
Gk(w) in (B3) is now replaced its finite-temperature generalization G\ (uj) (2.49). The outer region discussion 
is exactly the same as before except that one should use finite temperature counterparts of / and A t in (B4). 
For T <C /i, the outer region equation depends on T analytically and the asymptotic behavior in equations (B5) 
and (B6) is unchanged. Thus for T <C /i, one simply replaces Gk(w) in (Bll) by G^\uj), but keeping in mind 
that now a±,b± also have analytic temperature dependence. 



45 The Wronskian of equation (Bl) is given by 
which is independent of z. 



V^i,<£ 2 ] = -^(<l>id z <h -fadzfa) (B18) 
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Appendix C: Double trace deformations 



Consider a large TV gauge theory in Euclidean signature deformed by a double-trace operator 



SS 2 



X -k J d d xQ 2 (CI) 



where O is a scalar operator whose Euclidean correlation function (in the absence of deformation) is given by G(x). 
With this deformation the two point function for O now becomes: 

Gk = ^(o{x)O{0)e-^S dd y°^ 2 ), 

Z K = (e-^f dd y°M 2 ) (C2) 

which leads to 

G K (x) = ±- £ t££ (f[ [ d d y m ) (O(x)O(0)O( yi ) 2 . . . G(y n ) 2 ) (C3) 

^ n=0 Z U - \m=l J ) 

The disconnected diagrams cancel between the numerator and denominator, leaving us with only connected diagrams. 
To leading order in the large N limit, the n-th term in the above equation becomes 

{-n) n J d d Vl . . . d d y n G(x- yi )... G(y n -! - y n )G(y n ) . (C4) 

In momentum space we thus find a simple geometric sum 

GM = £<-«)-0*» W = Y f^ ¥ - ■ (05) 

n=0 x ' v ' 

Analytically continuing (C5) to the Lorentzian signature and using the standard relation between the Euclidean and 
retarded functions we conclude that 

G ( n\^k)= 1 (C6) 

G R (cj, k) + k 

where Gr(uj, k) and G^\lj, k) are retarded functions for O before and after the double trace deformation. 
For a CFT, the Euclidean correlator G(k) in momentum space for an operator O can be written as 

G(k) = C(A)k 2A ~ d (C7) 

where where A is the scaling dimension of O and C(A) is some constant factor. Plugging the above equation into (C5) 
we find that for A > |, 

G K {k)->G(k), k^O. (C8) 

This is expected, since for A > |, the double trace deformation (CI) is irrelevant and the theory flows back to the 
original fixed point in the IR. When A < |, for which (CI) is relevant, we find that in the IR 

G K (k) -^G-\k) + . . . , k -> . (C9) 

The first term in the above equation is a constant (which gives rise to a contact term in coordinate space) and can be 
ignored. The second term corresponds to an operator of dimension d — A. We thus conclude that under the double 
trace deformation the system flows to a new CFT in the IR in which O has dimension d — A. Note that two-point 
functions of other operators (which are orthogonal to O) are unaffected. 

In the bullet below that of equation (2.11) in Sec. (II A) we mentioned that for vjj G (0,1), there are two ways 
to quantize a scalar field <p in AdS. In the standard quantization (Dirichlet boundary condition) the corresponding 
boundary operator + has dimension A + = | + ujj, while in the alternative quantization (Neumann boundary 
condition), the corresponding boundary operator O- has dimension A_ = | — vjj. From our discussion above, upon 
turing on a double trace operator 2 _ in the alternative quantization, in the large N limit the system should flow in 
the IR to a new fixed point in which O- has dimension d — A_ = A + with the dimensions of all other operators 
unchanged. Clearly this new fixed point should be identified with that corresponding to the standard quantization. 
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Appendix D: Fourier transforms 



Consider an integral of the form 



I(x) = J dk F{y k y 



ikx 



(Dl) 



where v k = c^/k 2 + k$ with c and fco some constants. We will assume that the analytic structure of F is such that 
the only singularities in the complex /c-plane are at the beginning of the branch cuts in z/, that is at k = ±zfco- If 
that is so, then we can take the two branch cuts to go upwards and downwards in the complex /c-plane and then (for 
positive x) by deforming the contour upwards we can rewrite the above integrand as 



I(x) 



k 2 



k 2 

^0 



(D2) 



Now in the limit of large x we expect the integral to be dominated by the lower limit of the integrand n ~ fco- We 
can expand the sub-exponential part of the integrand in powers of (k — fco) and perform the integrals directly; thus 
it is clear that the leading x behavior will be of the form 



I(x) 



-kox 



(D3) 



where the power appearing in depends on the exact form of the integrand. The key point here is the the location 
of the branch cut sets a scale which determines the correlation length in the integral. 



Appendix E: Density of states of semi-local quantum liquids 



In this Appendix we consider the density of states of a semi-local quantum liquid. The discussion is motivated by 
a recent argument of [92], which appears to indicate that scaling symmetry only in the time direction cannot persist 
to arbitrarily low energies. The argument is very simple and goes as follows. On dimensional grounds, the scaling 
symmetry in the time direction implies that the density of states of such a theory should have the form 

p(E) = A5{E) + | (El) 

where A and B are constants. The term proportional to S(E) gives rise to an entropy density, while the term 
proportional to 1/E appears to be needed to have any nontrivial dynamics (as otherwise the Hamiltonian will be 
proportional to the identity operator in this sector). With ^ 0, however, the integral of the 1/E term is divergent 
at low energies, suggesting an infinite number of states at arbitrarily low energies. This peculiar result suggests that 
the scaling symmetry must always break down at sufficiently low energies, say below some new dynamically generated 
scale Eir. Note that the argument is independent of the existence of a ground state entropy density, i.e. it should 
apply even if A = 0. 46 

We note, however, there are some subtleties with the argument. In order to have a finite density of states we should 
introduce an IR cutoff, i.e. we should put the system in a finite box with volume Vd-i = L d ~ x , and (El) applies to 
the system in a finite box. Even without a new dynamically generated scale Em, in a finite box, we expect that the 
system should develop a finite gap e(L) with 

e(L) -> (E2) 

in the thermodynamic limit L — >• oo. Thus, more precisely, the argument for a new dynamically generated scale Em 
should be made in the thermodynamic limit, i.e. Em should remain finite for L — >• oo. Let us now assume that such 
an Em does not exist, and see whether one can find an immediate contradiction from (El) in the thermodynamic 
limit. 



For example, it will rule out the kind of local quantum critical point discussed in the literature (e.g. [29]) as a genuine low energy 
description of of certain quantum phase transitions in heavy fermion materials. 
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In general it is hard to deal with the total density of states of an interacting system. In the thermodynamic limit, 
we can instead consider entropy whose exponential can be considered as giving the total density of states. Thus we 
first translate equation (El) into a statement regarding the entropy. For E ^ 0, taking the logarithm of (El) we find 

S{E) =\ogB-\ogE . (E3) 

The \ogE term is singular in the E — >> limit, but note that this singular term is not proportional to the volume, and 
thus appears to be negligible in the thermodynamic limit. This by itself does not imply that (El) is fine with B ^ 0. 
But it does indicate that drawing a conclusion requires going beyond the thermodynamic limit, in which case one will 
need to worry about the precise boundary conditions, shape of the box, and etc. In particular, given that the finite 
volume cutoff breaks the scaling symmetry, it is not clear whether (El) is really valid. To further highlight possible 
subtleties associated with going outside the thermodynamic limit, note that there is an alternative way to obtain the 
entropy from (El), 

N(E) = / dE' p(E') = e s ^ E) . (E4) 
Jo 

In the thermodynamic limit this definition is equivalent to simply taking the logarithm of the density of states. But 
applying (E4) to (El) we find that 

N(E) = B log -+ S(£) = log5 + loglog E>e(L) (E5) 

where we have put an IR cutoff e(L) in the integration of (E4). The ^-dependence in (E5) differ from that in (E3). 
In the L — >• oo limit the second term in S(E) is singular, but the divergence is much smaller than volume factor in 
the usual thermodynamic limit. 

The above discussion does motivate a stronger argument based on the thermodynamic limit. Instead of (El) one 
may argue that the scaling symmetry in the time direction implies that the entropy of the system should have the 
form 

S(T) = V d . 1 (a + blogT) (E6) 

where for convenience we have used the canonical ensemble and a, b are T-independent constants. Argument similar 
to that below (El) again seems to say that nontrivial dynamics requires 6^0, which then implies that the system 
should have a dynamically generated cutoff in the thermodynamic limit. We note, however, the coefficient b of the 
logT term may not be required to be nonzero to have nontrivial dynamics, as the systems which exhibit such a scaling 
symmetry have a UV completion. The leading contribution could come from an irrelevant operator. We will discuss 
a potential example below. 

For a boundary theory at finite density described by a charged black hole in the bulk, the entropy has the form 

S(T) = V d -t (N 2 a(T) + b(T) + 0(1/7V 2 )) . (E7) 

The 0(N 2 ) term comes from the black hole geometry. The term proportional to b(T) comes from a one- loop calculation 
in the bulk. It is a well defined question to ask whether b(T) has a singular logT piece in the T — >• limit. Here we 
will not attempt a full one-loop calculation. To get an indication of the behavior of b(T), we compute instead the 
single-particle density of states^ 7 p s (uj) for a scalar field in AdS2 x Using p s one obtains an "off-shell" partition 

function log Z s (T) using the standard formulas 



n=l 



n 



oo 



log Z s {T) = YJ s [-\, f s (T) = I duj Ps (oj) (E8) 



where f 3 (T) is the single-particle partition function in the bulk. We would like to see whether there is a singular logT 
term in logZ s (T) which can be traced to a similar term in f s (T) 48 . 
Before going into technical details, let us summarize the results: 



It should be kept in mind this bulk single-particle density of states does not correspond to any single-particle density of states in the 
boundary. 

The standard one loop partition function at a finite T which can be used to directly extract b(T) is computed from path integral 
in the Euclidean black hole geometry with given T. There exists argument (see e.g. [147]) that the "off-shell" partition function we 
are considering is equivalent to that computed from Euclidean path integral. Clearly it would be desirable to compute the Euclidean 
partition function directly. 
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A. For a scalar in A0IS2 x introducing a finite volume in the field-theory spatial directions does not result in 
a finite density of states in the bulk 49 , due to the presence of the horizon. To proceed, we will thus introduce 
an additional cutoff e ~ A -1 where A is a radial coordinate cutoff of the spacetime above the horizon. 

B. Within this regularization scheme, for a scalar operator of real scaling dimension in SLQL, the dangerous 
log T piece is not present. Instead the leading cutoff-independent contribution in the density of states comes from 
an irrelevant operator (a double trace deformation) and is integrable at small uj. One finds that f 3 (T) ~ T 2vk =° 
and thus the corresponding contribution to logZ s (T) goes to zero in the T — >> limit. 

C. For an operator with complex scaling dimension in SLQL, one does find a non-vanishing (cutoff-independent) 1/uj 
piece in the single-particle density of states. Upon integration, this leads to a singular single-particle partition 
function f s ~ logT and thus a singular logZ s (T). We would conclude that in this case, a new dynamical scale 
Eir should emerge in the thermodynamic limit. This is in nice correspondence with the discussion of sec. HE 3, 
where we show that precisely in this case the system is unstable to the condensate of the scalar operator and 
thus a new dynamical scale must emerge. We expect similar results should also apply to a fermionic operator 
with a complex dimension, which will again be in nice correspondence with the fermionic instability discussed 
in Sec. HIE. 

From the above, it then appears that the argument based on (El) or (E6) alone is not enough to rule out the 
presence of a scaling symmetry in the time direction persisting to arbitrarily low energies; it seems that dynamics still 
plays a role in determining whether this is possible. Clearly further study of this is desirable. 



1. Single-particle density of states in the bulk 

Here we calculate the density of states of a neutral scalar field <fi in AdS2 x whose equation of motion is given 

by (A6)-(A8) with q* = 0. We will consider solutions which are normalizable at the boundary 50 , i.e. a = in (A9). 
At the horizon £ — >• 00, the correponding wave function is given by a plane wave 

</>(u>, fc, C) = e iujC+iS ^^ k) + e-^C-iSi^) ( E9 ) 

with 5((j,T/k) a phase shift. Note that the normalizable condition and a discrete k does not impose a quantization 
condition on uj; the system has a continuous spectrum (we are dealing with a scattering problem) even with a finite 
volume cutoff. While we expect that the system should develop a discrete spectrum in a finite volume, the spacings 
of low-lying spectrum could be much smaller than 1/L, with L is the system size, in the large A and large N limit, 
and thus not immediately visible at our level of approximation. Thus in order to obtain a finite density of states, we 
will impose an IR cutoff by hand at £ = A near the horizon 51 . A has units of time and is the inverse of the energy 
cutoff scale E c discussed above. This now requires that <\> satisfies 

0(cj,A) = O ujA + 5(u, 5 k ) = (n + tt, n = 0, 1, • • • (E10) 

Taking differential on both sides of the last equation we find that 

duo A+— )=7rdn -> — = - I A + — ) (Ell) 
\ duo J duo 7r \ duj J 

We thus conclude that the density of states is given by 



/nx, f dd ~ lk dn V d -i f d d - x k . p / \ /t~\ -t c\\ 

where f(oj) is the cutoff-independent part given by 

= —] (2^ ■ (E13) 

Now note that the scaling symmetry of AdS2 implies that the phase shift 5((j,Sk) is in fact uo independent as the only 
dependence of a wave function on uj should be through uj(. Thus (E13) is identically zero and we are left with a piece 
which depends on the cutoff. 



As emphasized earlier the relation of this bulk single-particle density of states with the density of states in the boundary theory is rather 

indirect. In particular, one should not view the cutoff e introduced below as a new emergent dynamical scale. As discussed in [147] the 

divergence of log Z s (T) in the e — > limit in fact is identical to UV divergences in the Euclidean partition functions. 

For defmiteness we will restrict to the standard quantization. 

A should be fc-dependent, but we will not specify the dependence. 
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a. Double trace deformation 

The above discussion is, however, not the full story as the AdS2 region appears in the IR part of the full black 
hole geometry, which provides the UV completion. As discussed in Sec. II D (see discussion around equation (2.60)), 
the leading irrelevant deformation of SLQL is provided by the double trace operator, which we now consider. With a 
double trace deformation, the normalizable solution is now given by setting 

a = £k(3 (E14) 

in (A9). Note that the double trace parameter has dimension —2^. Up to an overall scaling the normalizable 
condition can also be written as 

= ffcC*""* + C* +Vk C -> (E15) 
The full normalizable solution in AdS2 then has the form 

= &r(i - 1*) (2) v / CJ r -, fc K) + r(i + ^)(-J x/C^K) (Eie) 

We thus find the phase shift is given by 

1 + Ce-f**' U \2i r(l + i/ fc ) [ ' 

For small uo we thus find that 

5 (<u) = const + C sin u^tt + • • • (E18) 

which when plugged into (E13) leads to 

where • • • denotes higher order corrections in uj. The above density of states is integrable at uj = 0. 52 We thus find 
that that the cutoff-independent part of the density of states is controlled by leading irrelevant perturbation. 



2. Imaginary v 

Let us now consider an imaginary v = — iX. For notational simplicity we will now drop the momentum subscript. 
Embedding the AdS2 into the full extremal BH geometry) amounts to requiring the normalizable wave function to 
have the asymptotic behavior near the AdS2 boundary 

<j> = e - i9 C^ iX + e i9 C^- iX (E20) 

for some choice of (which depends on A), whose precise value will not concern us. The full normalizable solution in 
AdS2 is then almost identical to (El 6) with v now imaginary 

ab = e- ie T(l - v) (^Y^J_ u ((dz) + e i6 T(l + v) ^z~J u (ujz) 
= e- %e ^z~J lX (ujz) + e i@ ^J- iX ((dz) ~ (E21) 

with 

6 = A log (d + const (E22) 

Note that in contrast to the earlier case, is now dimensionless, which indicates that f((d) now has to be proportional 
to co> _1 on dimensional ground. The phase shift is now given by 

e 2 ^ = * (E23) 



Recall that for fe < /i, ^ ?s > 0- 
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We thus conclude that 

1 dS A e 2A7r - 1 



7r duo uj 1 + e 2X7r + e A7r cos 26 



(E24) 



which has an un-integrable singularity at uj = 0. Note that including the momentum integrations in (E13) will not 
cure this non-integrability as v\~ is almost k independent for k <C /i*. 



Appendix F: Review of critical exponents 



In the vicinity of a critical point we observe scaling behavior of various observable quantities, which is characterized 
by a set of critical exponents. We list some of the most commonly used exponents in the following. We will denote 
the external tuning parameter g with which we tune the system to the critical point g = g c . Near the critical point 
the spatial correlation length diverges as 

^-iff-ffcr^ • (fi) 

The energy gap for elementary excitations scales as 

Egap ~ C ~ \9 - 9cr ZVerit , (F2) 

where z is called the dynamic critical exponent and determines the characteristic time scale of the approach to 
equilibrium via r eq ~ 1/E gap . On the condensed side the order parameter O also exhibits scaling near the critical 
point; the corresponding exponent is: 

(0)~\g-g c f , (F3) 

and exactly at the critical point it will depend on the source as 

(O) ~ J 1 ' 8 , (F4) 

where the coupling to the external source is JO. The correlation function \ — (OO) can also be used to probe the 
physical properties of the critical point. The next critical exponent we introduce is for \ a ^ zero momentum: 

x(fc = , w = 0)~ \g-g c r . (F5) 

The decay of Gr at the critical point in the free theory would be l/x d ~ 2 , the deviation from this is characterized by 
rj: 

X ( x ,u, = 0)\ g=go ~-^. (F6) 
To study the scaling of thermodynamic functions we introduce a as: 



f~\g-g c \ 2 - a , (F7) 



where / is the free energy density. 
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